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ABSTRACT 

Analyses of spectra obtained with the Far Ultraviolet Spectroscopic Ex- 
plorer (FUSE) satellite, together with spectra from the Copernicus and IMAPS 
instruments, reveal an unexplained very wide range in the observed deu- 
terium/hydrogen (D/H) ratios for interstellar gas in the Galactic disk beyond 
the Local Bubble. We argue that spatial variations in the depletion of deuterium 
onto dust grains can explain these local variations in the observed gas-phase D/H 
ratios. We present a variable deuterium depletion model that naturally explains 
the constant measured values of D /H inside the Local Bubble, the wide range of 
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gas-phase D/H ratios observed in the intermediate regime (logiV(H I) = 19.2- 
20.7), and the low gas-phase D/H ratios observed at larger hydrogen column 
densities. We consider empirical tests of the deuterium depletion hypothesis: (i) 
correlations of gas-phase D/H ratios with depletions of the refractory metals iron 
and silicon, and (ii) correlation with the H2 rotational temperature. Both of 
these tests are consistent with deuterium depletion from the gas phase in cold, 
not recently shocked, regions of the ISM, and high gas-phase D/H ratios in gas 
that has been shocked or otherwise heated recently. We argue that the most 
representative value for the total (gas plus dust) D/H ratio within 1 kpc of the 
Sun is > 23.1 ±2.4 (1 a) parts per million (ppm). This ratio constrains Galactic 
chemical evolution models to have a very small deuterium astration factor, the 
ratio of primordial to total (D/H) ratio in the local region of the Galactic disk, 
which we estimate to be fa < 1.191q 15 (1°") or < 1-12 ± 0.14 (ler) depending on 
the adopted light element nuclear reaction rates. 

Subject headings: ISM: abundances — ISM: structure — Galaxy: abundances - 
Galaxy: solar neighborhood — Galaxy: deuterium — Ultraviolet: ISM 



1. INTRODUCTION 

Accurate measurements of the deuterium/hydrogen (D/H) ratio (by number) provide 
critically important tests for models of primordial nucleosynthesis, Galactic chemical evo- 
lution, and the chemical properties of the intergalactic medium (IGM). Among all of the 
light elements, deuterium arguably provides the most stringent tests of our understanding 
of these three critically important topics, because the evolution of deuterium is predicted to 
be rather simple and the dependence of D/H varies steeply with baryonic density. Epstein, 
Lattimer, & Schramm (1976) and others have argued that deuterium was formed only in 
the Big Bang and that D is easily converted to 3 He, 4 He, and heavier elements by nuclear 
reactions in stars, a process commonly called astration. Thus, as stars evolve, the interstellar 
medium (ISM) receives deuterium-depleted and metal-enhanced material from stellar winds 
and supernova explosions, leading to a monotonic decrease in the D/H ratio with increased 
processing through stars. Other light elements (e.g., 3 He, 4 He, 7 Li) have more complex 
creation and destruction histories or otherwise provide less sensitive tests of these theories 



1 Bascd on observations made with the NASA-CNES-CSA Far Ultraviolet Spectroscopic Explorer. FUSE 
is operated for NASA by the Johns Hopkins University under NASA contract NAS5-32985. 
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(cf. Boesgaard & Steigman 1985; Wilson & Rood 1994). Infall of less astrated gas, which 
is deuterium rich and metal poor, from the IGM and small galaxies captured by the Milky 
Way Galaxy, complicates models of Galactic chemical evolution, because the rate of infall, 
its time and spatial dependence, and the timescale for mixing with other gas are poorly 
known. 

Although based on different physical processes and relating to very different times in 
the early universe, the two main approaches for determining the primordial D/H ratio, 
(D/H) prim , are now in agreement. The first approach is to measure the column densities of 
D I and H I in quasar absorption line systems. Kirkman et al. (2003) have determined a 
mean value of D/H = 27.8^3'g parts per million (ppm) by averaging the log (D/H) values 
using column densities iV(H I) and iV(D I) obtained from absorption line measurements of 
gas in the lines of sight toward five quasars. Most of these measurements use high resolution 
spectra obtained with the Keck telescope. Since the observed gas has very low metallicity, 
these values of D/H are likely close to, but perhaps slightly smaller than, the primordial (Big 
Bang) value, (D/H) pr i m . The measurement of D/H = 22 ± 7 ppm for the high velocity cloud 
Complex C (Sembach et al. 2004) likely refers to gas located 6-10 kpc into the Galactic halo 
with low metallicity, 10-30% solar abundance (Collins et al. 2003), and this value should be 
somewhat smaller than (D/H) prim . 

The second approach involves the analysis of data from the Wilkinson Microwave 
Anisotropy Probe (WMAP) satellite and other cosmic microwave background experiments. 
Assuming a flat A-dominated universe and other standard assumptions, several authors 
(Spergel et al. 2003; Cyburt, Fields, & Olive 2003; Coc et al. 2004; Sanchez et al. 2006) have 
derived tight constraints on the baryon closure parameter, fl^, and the Hubble constant. 
Together with other constraints from the power spectrum of galaxy clustering and La forest 
data, they find that VL b h 2 = 0.0224 with a la error of about 4%. These results are summa- 
rized in Table 1. The ratio of baryons to photons (multiplied by 10 10 ), i]i , is simply related 
to Qt>h 2 (Buries, Nollett, & Turner 2001). In the standard big bang nucleosynthesis (SBBN) 
model, rjio is the single parameter that predicts the primordial light element abundances 
including (D/hl) prim . For example, Spergel et al. (2003), Coc et al. (2004), and Sanchez et 
al. (2006) obtain (D/H) prim rs 26 ppm. Using the same values of Q b h 2 and r] 10 but a more 
recent compilation of nuclear reaction rates (Angulo et al. 1999), Cyburt, Fields, & Olive 
(2003) obtain (D/H) prim = 27.5^;g ppm. The somewhat different values of (D/H) prim arise 
from the different authors treating the light element nuclear reaction rates and error propa- 
gation somewhat differently (cf. Nollett & Buries 2000). We therefore use in the subsequent 
analysis two different values for (D/H) prim , 27.5l^g ppm (Cyburt, Fields, & Olive 2003) and 
26.0+ 

L? PP m (Coc et al. 2004). Both values are consistent with the quasar absorption line 
result for (D/H) prim . 
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Beginning with observations by the Copernicus satellite, data from high-resolution ultra- 
violet spectrometers on several satellites have been analyzed to measure D/H in the Galactic 
ISM. These results provide insight into the abundance of deuterium in the present epoch, 
which with measurements of various metal abundances constrain models of Galactic chem- 
ical evolution. The International Ultraviolet Explorer (IUE) satellite and more recently the 
Goddard High Resolution Spectrometer ( GHRS) and Space Telescope Imaging Spectrograph 
(STIS) instruments on the Hubble Space Telescope (HST) have provided UV spectra of the 
Lya absorption lines of H I (1215.67 A) and D I (shifted by -82 km s" 1 relative to the H I 
line), which are needed to measure D/H in the ISM along lines of sight to stars located 
within a few hundred parsecs of the Sun. The same analysis technique has been applied to 
hydrogen and deuterium lines higher in the Lyman series that lie in the spectral range 912- 
1025 A, using data obtained with the Copernicus, the Interstellar Medium Absorption Profile 
Spectrograph (IMAPS), and the FUSE satellite. For a description of the FUSE satellite and 
its operations, see Moos et al. (2000) and Sahnow et al. (2000). 

The far-UV spectral range is essential for sampling more distant lines of sight using 
hot subdwarfs and white dwarfs as background sources, because the decreasing opacity in 
the higher Lyman lines allows one to observe the deuterium lines without masking by the 
nearby very opaque broad hydrogen lines. For example, the dark core of the H I Lya 
line absorption broadens to mask the D I Lya line when the hydrogen column density 
logiV(H I) > 18.7, but at Ly7 (972 A) the horizon for detecting the DI Ly7 line extends to 
logiV(H I) ps 20.0. Measurements of D/H using the higher Lyman series lines are feasible to 
about logiV(H I) ~ 21 (cf. Hoopes et al. 2003), where the density of overlapping H 2 lines 
and the large widths of the interstellar H I Lyman lines determine the final horizon for such 
analyses. 

There are several subtle effects that can enter the analysis of D I and H I lines in the far 
UV. One disadvantage is that, for long lines of sight, there is an increased chance of weak 
features at large velocity displacements. An H I feature near -82 km s" 1 can masquerade 
as D I absorption. An important advantage is that the far UV includes a large number of 
Lyman lines with a wide range of optical depths. Thus, measurements of iV(H I) and A^(D I) 
can be based on several spectral lines in the Lyman series. Another advantage is that, for the 
higher Lyman lines, the upper state lifetimes are long, decreasing the damping and better 
separating the H I and D I absorption for each line in the Lyman series and between adjacent 
lines in the series. 

Since the very first of these measurements, a disturbing aspect of the D/H studies has 
been the wide range of measured values with no apparent pattern in the data or a convincing 
explanation for it. We address the first part of this problem with the latest data, including 
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many measurements with FUSE] we conclude (see § 2) that the variations are real, and 
we show that there is a consistent pattern. We propose that the variations arise from time- 
dependent deuterium depletion onto dust grains (see § 4 and § 5), and we examine two critical 
tests (see § 6) of this model. The first test (§6.1) uses iron and silicon depletions as proxies 
for grain destruction, and the second test (§6.2) uses correlations of D/H with H2 excitation 
temperatures. We conclude (see § 7 and § 9) that the most likely value for the total D/H 
ratio within 1 kpc of the Sun is about 45% larger than the previously assumed value of «15 
ppm, which provides a revised constraint on models of Galactic chemical evolution. Also, 
the inferred depletion of deuterium from the gas phase places constraints on the physics and 
chemistry of interstellar dust. See Linsky (2006) for an early version of this paper. 



2. MEASUREMENTS OF D/H IN THE LOCAL GALACTIC DISK 

Using FUSE data, Wood et al. (2004) measured (D/H) gas along the lines of sight to two 
hot subdwarf stars, JL 9 and LS 1274, located at distances of 590 ± 160 pc and 580 ± 100 pc, 
respectively. They also summarized the (D/H) gas measurements observed by the previously 
mentioned satellites for 38 lines of sight to stars located between 3.2 pc and 2200 pc from the 
Sun. These column densities are summarized in Table 2 together with the new results for 
the PG 0038+199 (Williger et al. 2005), HD 90087 (Hebrard et al. 2005), LSE 44 (Friedman 
et al. 2006), WD 1034+001 and TD1 32709 (Oliveira et al. 2006) lines of sight and (D/H) gas 
for Lan 23 (Oliveira et al. 2003), which was not included in previous compilations owing 
to the large uncertainty in (D/H) gas . Table 3 lists the (D/H) gas ratios computed from the 
measured column densities in Table 2, and the la errors in (D/H) gas are computed from the 
square root of the sum of the squares of the percent errors in iV(D I) and N(R I). When 
authors present their results for (D/H) gas but not iV(D I), we have computed iV(D I) and its 
uncertainty to be consistent with the published values of (D/H) gas and iV(H I). The (D/H) gas 
measurements plotted in Figure 1 show a complex picture that, as we shall see, seriously 
challenges previously held ideas concerning Galactic chemical evolution. 

Table 2 lists the most recent measurements of iV(H I) and iV(D I), ordered by increasing 
iV(H I), for interstellar gas along the lines of sight to the listed stars, obtained using the 
FUSE, Copernicus, IMAPS, and GHRS and STIS instruments on HST. (D/H) gas , which is 
derived from the N(D I)/iV(H I) ratio, is not always the same as originally reported because 
there are measurements of iV(H I) obtained more recently with higher S/N and/or spectral 
resolution than the corresponding N(D I) measurements. This is typically the case when 
the only measurement of N(D I) was obtained with the Copernicus satellite and iV(H I) was 
later measured from HST spectra. For the sightlines to \x Col, (3 CMa, 9 Car, and A Sco, we 
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list iV(H I) and N(D I) for a specified velocity component (see notes to Table 2) rather than 
for the entire inhomogeneous line of sight, since the high quality spectra enabled the authors 
to identify the velocity components and measure the column densities in each component, 
although the components were likely not resolved. 

For the sightlines that extend through the Local Bubble to more distant stars, we 
also list (D/H) gas _LBj the value of (D/H) gas for the line of sight beyond the Local Bubble 
foreground, obtained by subtracting iV(H I)lb = 10 19 ' 20 cm~ 2 and N(D I)lb — 10 14 ' 39 cm~ 2 
from the measured values of iV(H I) and N(D I) for the full line of sight. We assume 
that the Local Bubble extends to logiV(H I) = 19.2, since beyond that point the (D/H) gas 
values are no longer consistent with a constant value. Sfeir et al. (1999) place the edge of 
the Local Bubble at logiV(H I) ~ 19.3, which corresponds to a 20 mA equivalent width of 
the Na I 5892 A line formed in cold interstellar gas. Within the Local Bubble, we adopt 
(D/H) gas = 15.6 ± 0.4 ppm (Wood et al. 2004). Given the complex three dimensional shape 
of the Local Bubble (Lallement et al. 2003), these assumptions for its column densities 
are approximations, but they provide good first estimates of the foreground required for 
isolating the deuterium and hydrogen column densities of the gas lying beyond the Local 
Bubble. This is especially important for lines of sight that do not extend very far beyond 
the Local Bubble, but this correction is small for the more distant lines of sight because 
the Local Bubble contribution to the total is minor. Those lines of sight that are definitely 
inside the Local Bubble are flagged as LBg in the eighth column of Table 2. 

Figure 1 shows the (D/H) gas ratios for the 47 lines of sight as functions of hydrogen 
column density. The (D/H) gas ratios extend over an unexpectedly wide range, but Wood et 
al. (2004) proposed that these results naturally fit into three regimes. For logiV(H I) < 19.2 
(in cm -2 units), corresponding to sightlines extending to the edge of the Local Bubble near 
100 pc, (D/H) gas has a constant value of 15.6 ± 0.4 ppm, where the quoted uncertainty is 
the la error in the mean for the 23 lines of sight within this regime. This constant value 
for the local (D/H) gas ratio had been previously established based on HST data (Linsky 
1998, 2003), but the addition of new lines of sight observed by FUSE (Moos et al. 2002) 
strongly supports this result. The very different values of (D/H) gas measured in the distant 
(logiV(H I) > 20.7) and intermediate (logiV(H I) = 19.2 — 20.7) regimes beyond the Local 
Bubble are described at the end of this section. 

Hebrard & Moos (2003), Hebrard et al. (2005) and others have derived (D/H) gas in 
an alternative way from the product (D/0) gas x (0/H) gas . Since the charge-exchange re- 
action between O I and H I is fast, the ionization of oxygen, hydrogen, and deuterium 
are very nearly the same. In the ionization equilibrium calculations for their best mod- 
els (Models 2 and 8) for the Local Interstellar Cloud, Frisch & Slavin (2003) find that 
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each element is about 30% ionized. Thus in the gas phase (O/H) = N(0 I)/iV(H I) = 
[N(0 T)+N(0 II)]/LV(H I)+iV(H II], and a similar relation can be written for (D/O). Em- 
pirically, the abundance of O relative to H shows little evidence for spatial variation (Meyer 
et al. 1998; Andre et al. 2003; Cartledge et al. 2004). Also, the outcomes for (D/0) gas 
may not suffer from the effects of line saturation, as there are O I lines that are optically 
thin. All three of these considerations should make Ola good proxy for H I. The (D/0) gas 
ratio should therefore behave in a very similar way to (D/H) gas . The problems with this 
technique are that uncertainties in two different measurements must be considered instead 
of one, and the amount of depletion of oxygen from the gas phase may depend on density 
(Cartledge et al. 2004). Table 4 summarizes the recent measurements of (0/H) gas . Oliveira 
et al. (2005) found that inside the Local Bubble, (0/H) gas is constant with a mean value 
of 345 ± 19 ppm. Meyer et al. (1998), Meyer (2001), and Andre et al. (2003) find slightly 
higher values of (0/H) gas for lines of sight extending beyond the Local Bubble with perhaps 
a slight decrease in (0/H) gas at large column densities (log N(H I)> 21.5). We assume 
(D/0) gas = (3.84 ± 0.16) x 10~ 2 (la), the mean value that Hebrard & Moos (2003) find for 
lines of sight inside the Local Bubble, and compute the values of (D/H) gas listed in Table 4. 
Inside the Local Bubble, using the value of (0/H) gas obtained by Oliveira et al. (2005), we 
compute (D/H) gas = 13.2 ± 0.9 ppm. For values of (0/H) gas derived from longer lines of 
sight (Andre et al. 2003), we find (D/H) gas = 13.9 ±0.7 to 15.7 ±0.9 ppm, depending on the 
number of lines of sight included. In their Models 2 and 8, Frisch & Slavin (2003) conclude 
that (0/H) gas = 380 ± 30 2 ppm, leading to (D/H) gas = 14.6 ± 1.3 ppm. These alternative 
ways of determining (D/H) gas are in agreement or smaller by 1-2<t compared with the direct 
measurement of (D/H) gas = 15.6 ± 0.4 ppm obtained by Wood et al. (2004). 

To put these numbers into perspective, we note the recent measurements of the solar 
ratio, (O/H) = 458 ± 53 ppm (Asplund et al. 2004) and 390 ± 63 ppm (Melendez 2004). 
Bahcall, Serenelli, & Basu (2005) argue that the Asplund et al. (2004) value for the solar 
oxygen abundance is too low to explain the helioseismological determinations of convection 
zone parameters and surface helium composition, but a recent increase in the estimated 
opacity near the base of the solar convective zone (Seaton & Badnell 2004) minimizes the 
problem even for the Melendez (2004) low abundance. In §3 we argue that about 120 ppm 
of oxygen could be depleted on to grains if Mg, Fe, and Si have solar abundances in the ISM, 
but Frisch & Slavin (2003) conclude that the ISM abundances of these elements may be only 
60-70% solar. We therefore take the measured solar value (458 ± 53 ppm) and that value less 
120 ppm as a plausible range in (0/H) gas . Then the plausible range in the inferred (D/H) gas 
= (D/0) gas x (O/H) is 13.0 ± 2.0 ppm to 17.6 ± 2.1 ppm. All of the Local Bubble values 



2 We arbitrarily assign an error of 8% to their value of (0/H) gas based on the discussion in their paper. 
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of (D/H) gas measured directly and indirectly (see Table 4) fall well within this range. The 
differences among these values and with the directly measured (D/H) gas = 15.6 ± 0.4 ppm 
may represent measurement uncertainty and perhaps some systematic effects. 

An important new result noted for (D/H) gas measurements by Hebrard & Moos (2003) 
and emphasized by Wood et al. (2004) is that sightlines with large N(R I) appear to have 
systematically low values of (D/H) gas . Wood et al. (2004) called attention to the clump 
of four low (D/H) gas values at logiV(H I) > 20.7, which we refer to as the distant regime. 
Whereas unambiguously low values for D/H in lines of sight extending beyond the Local 
Bubble have been known since the Copernicus results for 5 Ori, e Ori, and 9 Car (Jenkins 
et al. 1999; Laurent, Vidal-Madjar, & York 1979; Allen, Jenkins, & Snow 1992), the new 
results for JL 9 and LSS 1274, together with the (D/H) gas values determined from FUSE 
observations of HD 195965 (794 ± 200 pc) and HD 191877 (2200 ± 550 pc) from Hoopes et 
al. (2003) and HD 90087 (2740 ± 800 pc) from Hebrard et al. (2005) indicate an emerging 
pattern. For these five long lines of sight, (D/H) gas = 8.6 ± 0.8 (standard deviation of the 
mean) ppm, nearly a factor of two below the Local Bubble value. 

The intermediate regime for log N(H I) = 19.2-20.7 exhibits a factor of 4-5 range in 
(D/H) gas : lines of sight to five stars have high values (near 22 ppm), five stars have low values 
(near 6 ppm), and nine stars have intermediate values. Since (i) the Copernicus, IMAPS 
and FUSE data all show this wide spread in (D/H) gas values in the intermediate regime, (ii) 
there does not appear to be any correlation of (D/H) gas with Galactic latitude, longitude, 
or distance, and (iii) the (D/0) gas values also show a large spread, we conclude that the 
wide spread in (D/H) gas in the intermediate regime is not due to instrumental artifacts or 
unknown systematic errors, but is rather a robust result that requires a sensible scientific 
explanation. We call attention to the five high data points in the intermediate regime with 
(D/H) gas PS 22 that will play an important role in our deuterium depletion model (see §5). 

3. PREVIOUS EXPLANATIONS OF THE D/H DATA 

The constant value of (D/H) gas inside the Local Bubble suggests that this matter has 
experienced a similar history of nuclear processing, mixing, and deuterium depletion (cf. 
§ 5). This history consists of star formation in the Scorpio-Centaurus Association some 
10 Myr ago and the creation of the Local Bubble by supernovae in the Lower Centau- 
rus Crux subgroup of the Scorpio-Centaurus Association (Maiz-Apellaniz 2001; Berghofer 
& Breitschwerdt 2002). Outside the Local Bubble, interstellar gas could have experienced 
different chemical evolution resulting from different rates of astration, infall of nearly pri- 
mordial material from the Galactic halo or the IGM, and insufficient mixing to homogenize 



- 9- 



the Galactic ISM. However, Galactic chemical evolution models are faced with the severe 
challenge of explaining the factor of 4-5 range in (D/H) gas ratios just outside of the Local 
Bubble and the clumping of low (D/H) gas values for the lines of sight with the largest val- 
ues of iV(H I). Recent Galactic chemical evolution models constructed to fit a wide range 
of constraints, including the abundances of deuterium, 3 He, 4 He, Li, and heavier elements 
created by astration, generally predict a moderate astration factor for the total amount of 
(D/H) tota i in all forms, i d = (D/H) prim /(D/H) tota j 1.5 (cf. Chiappini, Renda, & Matteucci 
2002; Romano et al. 2003). If we assume (D/H) prim = 27.5 ppm, then the astration factors 
in the intermediate regime mentioned in the previous section would lie in the range of 1.25 
to 4.6. How could such a large range in astration factors occur over distance scales of only 
100-500 pc, but an apparently constant astration factor of 3.2 appear to be representative 
for longer distance scales? 

Galactic chemical evolution models generally include enrichment of gas in the Galactic 
disk by the infall of less-astrated material from the IGM or gas from neighboring small 
galaxies. For example, Geiss, Gloeckler, & Charbonnel (2002) find that the Local Bubble 
value of (D/H) gas could be explained by a mixture of 60% gas representative of the protosolar 
nebula (with additional chemical processing over the last 4.6 Gyr) and 40% infalling gas 
with high (D/H) gas and low metals representative of the Large Magellanic Cloud. While 
this mechanism may explain the Local Bubble abundances, variations of the percentages to 
explain the wide range of observed (D/H) gas values would require that the infalling gas be 
spatially concentrated into small regions of the Galactic disk and then not be mixed. We 
note here that typical sizes for structures in the nearby infalling high-velocity cloud Complex 
C are about 100 pc (Collins, Shull, & Giroux 2003). We discuss the infall hypothesis further 
in §8. 

Both spatially variable astration and infall models may, in principle, explain the (D/H) gas 
measurements, but they require extreme assumptions concerning the deuterium destruc- 
tion and infall rates to explain the observed (D/H) gas ratios. Simulations with their three- 
dimensional supernova-driven ISM model led de Avillez & Mac Low (2002) to conclude that 
for the Galactic supernova rate, it takes some 350 Myr to completely mix the ISM on distance 
scales from parsecs to kiloparsecs. In view of this relatively rapid mixing rate, we expect 
the total interstellar D/H ratio in all forms, (D/H) to tai, to show little variation on length 
scales less than 1 kpc, unless one invokes localized starbursts sufficiently strong for astration 
to produce significantly lower (D/H) to tai on a timescale shorter than 350 Myr. Both the 
variable astration and infall models fail to explain the absence of any correlation between 
the observed (D/H) gas and (0/H) gas ratios (e.g., Draine 2005), which should show a negative 
correlation for both models. We explore instead another physical process that could have a 
much larger effect on local values of the (D/H) gas ratio than either astration or infall. 
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Figure 2 shows observed values of (D/H) gas and (0/H) gas for all sightlines that have 
published column densities of H, D, and O with la uncertainties of less than 0.1 dex for each. 
The error bars in the figure are la; because N(H I + 2H 2 ) appears in the denominator of 
both abscissa and ordinate, errors in (O/H) and (D/H) are partially correlated. Variations in 
astration or infall would produce anticorrelation of (D/H) total an d (O/H) total- Figure 2 shows 
the expected track of (D/H) to tai vs. (0/H) gas for chemical evolution model C-I of Chiappini 
et al. (2002), where (i) we have taken Chiappini et al.'s variation with galactocentric radius 
as a proxy for variations in astration and infall, (ii) we have rescaled the predicted (D/H) tota i 
values slightly to be consistent with primordial (D/H) to tai ~27.5 ppm estimated from quasar 
absorption lines and WMAP, and (iii) we assume partial depletion of O into grains (see 
below), with (0/H) gas = 0.75(O/H) to tai- The observed (D/H) gas values fall systematically 
below the (D/H) tota i values expected for at least this astration/infall model. The full (factor 
of ~3) range of variation in (D/H) gas is seen among the 6 sightlines with (0/H) gas between 
300 and 390 ppm. Variable astration and infall may account for a fraction of the variation 
in (D/H) gas , but it appears that most of the variation must be due to some other mechanism 
— such as depletion of D onto dust grains. 

The median value of (0/H) gas in Figure 2 is 363 ppm, but it is likely that another 
~ 120 ppm is locked up in grains in silicates or oxides (assuming the grains contain ~90% 
of solar abundances of Mg, Fe, and Si, with composition Mga;Fe2-iESi04). 3 From these well- 
characterized sightlines, we therefore estimate (0/H) total 483 ppm in the local ISM, slightly 
higher than but consistent with the recent estimate for the solar abundance (O/H) = 
458±53 ppm (Asplund et al. 2004). Thus we estimate that O gas /O tota i « (363/483) = 0.75 for 
diffuse regions. Note that for a primordial D/H = 27 ppm and (0/H) gas /(0/H) to tai ~ 0.75, 
the Chiappini et al. astration + infall model C-I (see Figure 2) has (D/H) to tai reduced to 
23 ppm by the time (0/H) gas has risen to ~ 360 ppm, very close to the measured value of 
345 ± 19 ppm (Oliveira et al. 2005) for the solar neighborhood. We will return to this point 
in §8. 



4. DEUTERIUM DEPLETION ONTO DUST GRAINS 

It is important to recognize that all of the direct (D/H) gas measurements refer to inter- 
stellar gas, whereas the Galactic chemical evolution and primordial nucleosynthesis models 
refer to (D/H) in all forms, (D/H) to t a i- Dust grains exist in both the cold (10-10 2 K) and 
warm (10 3 -10 4 K) phases of the ISM. In diffuse interstellar clouds, 50-90% of elements 



3 The abundance of H2O is negligible outside of dark clouds. 
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such as Mg, Si, and Fe reside in grains (e.g., Sembach & Savage 1996; Jenkins 2004), and 
gas-phase abundances of such elements relative to hydrogen are much lower than the total 
elemental abundances. A similar process may preferentially deplete deuterium in the gas 
phase relative to hydrogen. This could lead to large spatial variations in the gas-phase ratio, 
(D/H) gas , while the total deuterium abundance, (D/H) tota i, including both gas and dust, re- 
mains roughly constant. If deuterium depletion onto dust grains is important, then the total 
Galactic D/H ratio cannot be extracted from gas-phase measurements that ignore deuterium 
locked in dust grains. 

Jura (1982) first proposed that the depletion of deuterium onto dust grains might explain 
the factor of two difference in (D/H) gas values for the lines of sight to ( Pup and S Ori obtained 
from the analysis of Copernicus spectra (Vidal-Madjar et al. 1977; Laurent, Vidal-Madjar, 
& York 1979). Jura argued that the sticking probability of deuterium onto dust grains 
should be high and that deuterium should be more tightly bound to the dust than hydrogen, 
since the zero-point energies of deuterium-metal bonds are lower than for the corresponding 
hydrogen- metal bonds. Numerical calculations of the chemistry of grain mantles by Tielens 
(1983) showed that deuterium can be highly enriched in grains, because the lower zero- 
point energy for deuterium bonds as compared to corresponding hydrogen bonds strongly 
favors the formation of deuterated molecules like HDCO on grain surfaces (cf. Turner 1990), 
effectively removing deuterium from the gas phase in dark regions where grains are coated 
with ices. 

Draine (2004, 2006) developed this idea further by noting that the C-D bond energy 
in polycyclic aromatic hydrocarbon (PAH) molecules is 0.083 eV larger than the C-H bond. 
If H and D in the gas were in thermodynamic equilibrium with deuterated PAH material, 
one would have (D/H) dust /(D/H) gas « e 970 ^/ Td -% which exceeds 5 x 10 4 for T dust < 90 K. 
Although thermodynamic equilibrium arguments must be regarded with suspicion in the 
nonequilibrium ISM, Draine argued that it is plausible to consider that in a steady-state the 
grains might have (D/H) dust /(D/H) gas as high as 5 x 10 4 or more. 

Contemporary grain models require a substantial amount of hydrogen-bearing carbona- 
ceous material. From the relative strengths and profiles of different C-H stretching modes, 
Pendleton & Allamandola (2002) estimate that ~85% of the carbon in dust is in aro- 
matic form, with H/C~0.35, and ~15% of the carbon is aliphatic, with H/C~2.1. If (C 
in grains) /(H total) m 200 ppm, as estimated by Draine (2004), then the aromatic hy- 
drocarbon material would contain (H in aromatics) / (H total) ~ 60 ppm. Suppose that in 
the local Galactic disk (D/H) tota j = 22 ppm, and that in the aromatic material the D/H 
ratio is enhanced so that (D/H) dust = 0.27. Then (D/H) gas = (D/H) total - (D/H) dust x 
(H) dust /Htotai ~ 6 ppm. Therefore we see that if D enrichment can bring the aromatic grains 
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to (D/H)dust ~ 0.27, the gas phase deuterium abundance could be reduced to ~ 6 ppm, 
which would be sufficient to explain the lowest values of (D/H) gas observed. This would re- 
quire extreme enrichment of the grains, relative to the gas, by a factor (D/H)d us t/ (D/H) gas m 
0.27/6 ppm rs 5 x 10 4 . This is possible for Td ust < 90 K (see above). If the aliphatic and aro- 
matic hydrocarbons are equally D-enriched, then (D/H)d us t/(D/H) gas ~ 2 x 10 4 would suffice 
to explain the sightlines with the lowest values of (D/H) gas 6 ppm, with (D/H)d us t ~ 0.13. 

Figure 1 shows that (D/H) gas varies from one sightline to another, covering a range from 
~6 ppm to ~22 ppm. Draine (2004, 2006) and Wood et al. (2004) suggest a dynamic model 
of the ISM in which deuterium is depleted from the gas phase onto dust grains over time until 
the D-bearing grains are eroded or destroyed by sputtering or grain-grain collisions in strong 
shocks from supernova remnants, or by strong UV radiation fields from nearby hot stars, 
thereby returning the sequestered deuterium back to the gas phase. Is this idea supported 
by observations, and can it explain the three regime interpretation for the (D/H) gas data 
shown in Fig. 1? 

The most direct test of the deuterium depletion model would be to measure (D/H) dus t i n 
interstellar dust grains. Keller, Messenger, & Bradley (2000) have performed this critical ex- 
periment by measuring the deuterium abundance in interplanetary dust particles containing 
amorphous carbonaceous material. They argue that sample L2009*E2, which was captured 
in the upper atmosphere of the Earth, originated in the ISM because of the very high D/H 
ratios measured for various molecules (which indicates formation at low temperatures where 
fractionation is important), the elevated 15 N/ 14 N ratio, and the particle size distribution. 
For the whole dust grain sample, they find that (D/H)d U st = 220 ppm, but for one fragment 
of the sample (D/H) dus t = 1900 ppm. If we assume that the local ISM value of (D/H) gas is 
15 ppm, then the (D/H) dus t enhancement factors are 15 for the entire sample and 130 for 
the fragment. This measurement alone provides a proof of concept that deuterium depletion 
onto grains does occur in the ISM, although the highest D/H ratio measured by Keller, 
Messenger, & Bradley (2000) in this particle sample falls a factor of ~ 70 short of the factor 
of ~ 10 4 that is required for grain deuteration to explain the typical observed variations in 
(D/H) gas . It is interesting that high values of (D/H)d us t are only detected in amorphous or- 
ganic (i.e., carbonaceous) material in extraterrestrial samples (Keller, Messenger, & Bradley 
2000; Peeters et al. 2004). 

As part of their paper summarizing the observational evidence for deuterated PAHs 
(hereafter called PADs) in the ISM, Peeters et al. (2004) point out that deuterium fraction- 
ations (the ratio of deuterated molecules to undeuterated molecules) for simple molecular 
species are typically 0.01 to 0.1 (e.g., Markwick, Charnley & Millar 2001) and higher deu- 
terium fractionations have been measured (e.g., Parise et al. 2002). Even triply deuterated 
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ammonia (van der Tak et al. 2002) and methanol (Parise et al. 2004) have even been ob- 
served in cold interstellar clouds. Peeters et al. (2004) then argue that PADs are more stable, 
more abundant in the ISM compared to many of the simpler molecules that show deuterium 
enrichments, and carry more H atoms per molecule. PADs may, therefore, represent a large 
reservoir of deuterium-enriched species. 

Peeters et al. (2004) then call attention to infrared emission features at 4.4 /im and 
4.65 /im in ISO spectra of the Orion Bar and M17 and identify these features as C-D 
stretching modes in PADs. From the relative strengths of the emission bands of PADs and 
PAHs, Peeters et al. (2004) estimate that the D/H ratio in these molecules is 0.17 ± 0.03 
for the Orion Bar and 0.36 ± 0.08 for M17. These ratios are about 10 4 times larger than for 
atomic deuterium and hydrogen in the gas phase, (D/H) gas . It is striking that the D/H ratios 
required to account for the observed 4.4 /im and 4.65 /im emission are comparable to the 
degree of deuteration of PAH material argued for by Draine (2004, 2006), and which imply 
substantial depletion of D from the gas if this D/H ratio were characteristic of all of the 
hydrocarbon material in grains. Thus, deuteration of PAHs over time in a cold environment 
could be an important process for removing deuterium from the gas phase, and Peeters et 
al. (2004) describe four different astrochemical processes by which this could occur. 

Grains must be cold for the deuterium depletion process to operate on their surfaces. 
Even in the 7000 K gas in warm interstellar clouds inside the Local Bubble, the grains are 
cold because at the densities of diffuse clouds, grain heating is by absorption of UV photons 
from stars rather than by collisions with the gas; the grains cool by the emission of infrared 
radiation. The observed IR emission spectrum for wavelengths greater than 80 microns 
requires that typical grain temperatures are ~18 K (Draine 2003). 

Because (D/H) gas varies from one sightline to another, D depletion and return of D to 
the gas must be dynamic processes. What is the timescale for D depletion? Draine (2006) 
considers reaction of D + ions with neutral PAHs, and reaction of D atoms with PAH + 
ions, arguing that the impinging D will become chemically bound in the PAH, displacing a 
previously-bound H if necessary. Considering the charge distribution of the PAH population 
under "cold neutral medium" conditions (n H ~ 30 cm" 3 , n e 0.03 cm -3 , T ps 10 2 K), 
Draine (2006) estimates that the timescale for a D atom to collide with a charged PAH is 
~2 Myr. If no processes act to return D to the gas phase, D depletion by a factor of e could 
therefore be achieved in only 2 Myr. In more diffuse regions, the timescale will be longer. 
For "warm neutral medium" conditions («h ~ 0.3cm -3 , n e ~ 0.03cm" 3 , T « 6000 K), a 
similar calculation gives a much longer D depletion time of ~50 Myr. Therefore it appears 
that D depletion will be dominated by the cooler, denser regions where gas-grain collision 
rates are high. 
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An important point is that most of the surface area available for deuterium depletion 
is provided by the population of free-flying PAHs, but this population is thought to contain 
only CpAH/-f^totai ~ 60 ppm (Li & Draine 2001), or ~l/3 of the total aromatic material, with 
the remaining carbonaceous material in grains with radii a > 150 A. How can the material 
in the larger grains become D-enriched? In denser regions, free-flying PAHs are expected 
to coagulate and become part of larger grains. The population of small PAHs is presumed 
to be replenished by grain fragmentation in grain-grain collisions. Thus there is steady 
interchange of material between PAHs and large grains. The timescale for a PAH to collide 
and coagulate with larger grains is ~10 Myr in the cold neutral medium, assuming that the 
larger grains have a total projected area per H nucleon of 10 _21 cm 2 , and a characteristic 
velocity relative to the gas of ~ 1 km/s, driven by MHD turbulence (Yan et al. 2004). In 
warm neutral medium conditions, this timescale becomes very long; once again, it appears 
that the exchange of matter between PAHs and larger grains is dominated by processes in 
the cold neutral medium. Therefore, the conditions are favorable for the time-dependent 
depletion of deuterium from the gas phase due to local events without changing (D/H) to tai 
significantly by either astration or infall of less-astrated gas. 



5. TIME-DEPENDENT DEUTERIUM DEPLETION MODEL 

We now apply the deuterium depletion model to the three-regime interpretation of the 
(D/H) gas data presented in Figure 1. The measurements of (D/H) gas in the Local Bubble are 
for gas in warm clouds (T 7, 000 K) embedded in a hot gas bubble created by supernova 
events, which likely occured in the Lower Centaurus Crux subgroup of the Scorpio-Centaurus 
Association during the last 10 Myr (Maiz-Apellaniz 2001). The Local Bubble was last 
reheated and presumably shocked by a supernova event 1-2 Myr ago (Lyu & Bruhweiler 
1996; Berghofer & Breitschwerdt 2002). If the Local Interstellar Cloud is representative 
(Redfield & Linsky 2000), typical column densities for the warm clouds in the Local Bubble 
are iV(H I) ~ 10 18 cm -2 (see Table 2). The origin of these clouds is uncertain. It is possible 
that they consist of material that has cooled after being heated by thermal conduction or 
turbulent mixing with hot gas. In this case some fraction of deuterium in grains might have 
been returned to the gas by thermal sputtering while the gas was sufficiently hot for this to be 
effective. Without knowledge of the thermal history of individual Local Bubble gas clouds, it 
is not possible to predict what values of (D/H) gas are expected in this interpretation, but the 
observed value of ~15 ppm appears to be compatible with this picture. Alternatively, pre- 
existing neutral clouds might have been shocked by low- velocity shocks driven into the clouds 
by a sudden increase in ambient pressure. Frisch et al. (1999) argue that the abundances of 
Si, Mg, and Fe in the Local Interstellar Cloud are consistent with grain processing through a 
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shock with ~ 100 km s -1 . In this interpretation, grain-grain collisions would be expected to 
result in destruction of only a few percent of the refractory grain material like iron (Frisch & 
Slavin 2003). Since the condensation temperature (Savage & Sembach 1996) of Fe (1300 K) 
is much larger than for C (~ 50 K) to which the D is bound, we anticipate that grain-grain 
collisions, radiation, or other heating processes would strip D from grains before Fe. Since, 
as we will later argue, the total (gas and dust) value of (D/H) in the local Galactic disk is 
about 45% larger than (D/H) L B g asi we surmize that the shocks in the Local Bubble have only 
partially evaporated the deuterium-bearing grains that could have been formed in molecular 
clouds before the onset of star formation in the Association. As we shall see in §6.1, the 
modest depletion of deuterium in the Local Bubble corresponds to modest depletions of the 
refractory metals iron and silicon and, as described in §2, also oxygen. 

Radio observations in the H I 21 cm line provide the best information on the properties 
and number density of individual diffuse H I clouds in the ISM. Typical diffuse H I clouds 
have column densities in the range logiV(H I) = 18.8-19.8 with a median column density 
log iV(H I) = 19.3 (Lockman, Jahoda, & McCammon 1986; Lockman 2002). Cox & Reynolds 
(1987) cite a mean temperature for diffuse H I clouds of ~80 K. They also say that, on 
average, the number of such clouds per kpc in the disk with column densities in excess of 
iV(H I) is 5.7[iV(H I) / 10 20 cm" 2 ] -0,8 . Using this formula, the number of clouds per kpc with 
logiV(H I) > 19.3 is about 20, with a typical cloud separation of about 50 pc; however, the 
distribution should be very patchy. 

Lines of sight extending beyond the Local Bubble should pass though either recently 
shocked gas, not recently shocked gas, or a mixture of such regions. According to our 
time-dependent depletion model, the highest measured values of (D/H) gas should indicate 
recently shocked regions where (D/H) gas < (D/H) total , and the lines of sight with low values 
of (D/H) gas should indicate regions that have not been shocked for a long time and thus 
have severe deuterium depletion. Lines of sight extending beyond the Local Bubble with the 
highest measured values of (D/H) gas _LB, that is values above 21 ppm, include the lines of 
sight toward 7 2 Vel, Lan 23, WD 1034+001, Feige 110, and LSE 44. 

Given the patchy distribution of diffuse H I clouds, lines of sight extending beyond the 
Local Bubble can include none, one, or many of these clouds. If there is, at most, one diffuse 
H I cloud in a line of sight, then (D/H) gas should reflect the conditions of warm interstellar 
gas near the Sun that is partially ionized and likely has not had sufficient time to deplete 
an appreciable fraction of deuterium onto dust after the last shock. However, when a line 
of sight passes through several or many diffuse H I clouds, the composite value of (D/H) gas 
should be dominated by the abundances in these clouds where gas-phase deuterium can be 
highly depleted onto the cold grains even when the gas is not cold. Thus lines of sight in the 
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intermediate regime (logiV(H I) = 19.2-20.7) could contain little or no unshocked H I gas 
or be dominated by such material, leading to a wide range of (D/H) gas values. 

HST and Copernicus spectra of 7 2 Vel and diffuse Ha emission near the star show highly 
ionized gas (including O VI) along its line of sight from which Fitzpatrick & Spitzer (1994) 
and Reynoso & Dubner (1997) infer an H II region of path length 100-150 pc produced 
by a wind-blown superbubble and photoionizing radiation from 7 2 Vel and other massive 
stars. These conditions of strong ultraviolet radiation and recent supernova-produced shocks 
could remove deuterium from the mantles of dust grains, leading to the observed high value 
of (D/H) gas . The ionization along the other lines of sight with high (D/H) gas is not yet 
characterized. 

For lines of sight with larger column densities (logiV(H I) > 20.7), we expect that most 
of the column density will be contributed by cool H I gas, and the overall (D/H) gas value 
should be close to the low (D/H) gas values prevalent in the cool H I phase. Measurements of 
(D/H) gas along as many long lines of sight as possible with FUSE and other instruments are 
needed to test this hypothesis. For a discussion of possible line saturation effects see § 6.3.2. 



6. TESTING THE DEUTERIUM DEPLETION MODEL 

We now pose two tests to determine whether or not deuterium depletion is a valid model 
for explaining the wide range of (D/H) gas values for lines of sight extending beyond the Local 
Bubble. While neither test is conclusive, the statistical trends identified by these tests taken 
together make for a strong case. One problem in this analysis is that the lines of sight to the 
very bright stars analyzed with data from Copernicus have larger uncertainties in (D/H) gas 
than the lines of sight to a different group of stars studied using IMAPS or FUSE spectra. 
IMAPS reobserved the Copernicus stars 7 2 Vel, ( Pup, and 5 Ori, but FUSE was unable 
to reobserve any of the Copernicus stars to obtain more accurate D/H measurements. A 
rough estimate of the improvement in the measurements obtained with IMAPS and FUSE 
comes from a comparison of the (D/H) gas values for 7 2 Vel: the Copernicus result (York & 
Rogerson 1976) was 17.8lg° 6 4 ppm, whereas the IMAPS result (Sonneborn et al. 2000), which 
we use in our analysis, was 21.8±2.1 ppm. Although the two results are consistent within the 
Copernicus error bars, the IMAPS value is 2.1er (IMAPS error bars) above the value obtained 
from Copernicus spectra. Vidal-Madjar et al. (1977) used Copernicus data to infer two 
velocity components in the line of sight to ( Pup: in component 1 (D/H) gas = 22 ppm with 
no error estimate, and in component 2 (D/H) gas = 23^ ppm. The high spectral resolution 
and S/N of the IMAPS data allowed Sonneborn et al. (2000) to resolve or nearly resolve the 
complex velocity structure for this line of sight and to infer (D/H) gas = 14.0±2.3 ppm, which 
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is well below the Copernicus result. On the other hand, for the 5 Ori line of sight there is 
good agreement between the Copernicus result of 7.5 ± 2.3 ppm (Laurent, Vidal-Madjar, & 
York 1979) and the IMAPS result of 7.4±£J ppm (Jenkins et al. 1999). These comparisons 
suggest that we should be careful when including Copernicus data together with the more 
accurate recent data. In particular, the Copernicus data for the lines of sight to a Vir and 
ji Col have uncertainties in (D/H) gas that are near 100%. We have chosen to include all of 
the Copernicus data, which will clutter the subsequent figures somewhat, but in computing 
correlations we weight each of the data points by the inverse of the 1 a range in (D/H) gas . In 
this way, the recent data with small measurement uncertainties will dominate the correlation 
fits. 

In Figures 3-8 we compare (D/H) gas vs. other observables along the same lines of sight. 
Except for Figure 5, we do not use (D/H) gas _ L B because we have no reliable way of removing 
the Local Bubble foreground from the other quantity that (D/H) gas is being compared with 
for most lines of sight. For all of the correlation tests, we compute weighted least-squares 
linear fits, log[iV(D I)/iV(H I)] = a + bD(Fe) or a + bD(Si), to the observed data points 
by allowing all of the data quantities, JV(H I), N(D I), and either iV(Fe II) or JV(Si II), 
to vary randomly within their la error range. We compute a and b parameters for 1000 
such realizations and then fit Gaussians to the distributions of the a and b parameters. 
Table 5 lists the mean values of the a and b parameters and their la uncertainties, which 
are the half-widths of the corresponding Gaussian distributions. The relatively small la 
uncertainties in the b slope parameters for these fits demonstrates that correlations are real 
despite the correlation of the errors in the x-axis and y-axis quantities (see §6.3.6 for a 
different approach). For each plot we also apply the Spearman rank correlation test 4 , which 
provides an unbiased means for determining whether the null hypothesis of no correlation can 
be rejected for the unweighted data points. Column 9 in Table 5 gives the Spearman rank 
correlation coefficient, and column 10 gives the two-sided significance of its deviation from 
zero. Small values of this significance parameter indicate significant correlation. The quantity 
in parenthesis is the number of standard deviations by which the significance parameter 
deviates from its null-hypothesis value. These quantities are listed in Table 5 for the various 
tests described below. 



4 We use the IDL library routine FLCORRELATE to compute the Spearman (p) rank correlation test. 
The description of this test is from the IDL Reference Guide. 



-18- 



6.1. First Test: Correlation with Refractory Metal Depletions 

Draine (2004, 2006) proposed that a good test of whether deuterium depletion onto dust 
is responsible for the observed variations in (D/H) gas is to see whether there is a positive 
correlation of depletion of metals like Fe, Si, and Ti with (D/H) gas . Therefore, our first 
test of the validity of the deuterium depletion model is to determine whether for the same 
lines of sight, the observed (D/H) gas ratios correlate with the depletion of refractory metals, 
that is, elements that easily condense onto interstellar grains because of their high chemical 
bonding energies. In the low density ISM, the primary source of grain erosion is due to high- 
energy-sputtering processes that occur in supernova-generated shock waves (Jones 2000). 
Since iron in grains is very resistant to sputtering, small values of the iron depletion would 
indicate significant processing of the grains that, according to the model, would predict that 
most or all of the deuterium has been removed from the grains. Gas-phase iron abundances 
can be measured accurately in the ISM, as STIS and FUSE spectra both provide accurate 
measurements of the dominant ionization stage, Fe + , for both cold and warm gas. The 
depletion of iron, D(Fe) = log 10 [(Fe/H) gas /(Fe/H) Q ], varies from typically —1.2 in warm 
diffuse clouds to typically —2.2 in cool diffuse clouds (Savage & Sembach 1996; Jenkins 
2004). Very small iron depletions would indicate that the dust grains have been vaporized, 
depositing iron and deuterium from the dust back into the gas phase. Thus, very large 
negative values of D(Fe) should correlate with low (D/H) gas , if the deuterium depletion 
model is valid. Conversely, small negative values of D(Fe) should correlate with high (D/H) gas 
measurements. Column (8) in Table 3 lists the values of D(Fe) obtained using iV(Fe II) either 
from the reference given in column (11) of Table 2 or from the data compilations of Jenkins, 
Savage, & Spitzer (1986) and Redfield & Linsky (2002). In a few cases, we have revised 
slightly the values of iV(Fe II) as a result of more accurately determined f-values (Howk et 
al. 2000) for the optically thin lines that most strongly constrain the Fe II column density. 
These revisions range from +0.04 dex for BD+39°3226 to -0.04 for 9 Car. We then compute 
D(Fe) using our tabulated value of iV(H I). The undepleted value for iron is assumed to be 
the solar abundance ratio, log (Fe/H) = —4.55 ± 0.05 (Asplund et al. 2005), determined 
using a 3D hydrodynamical model of the solar atmosphere. 

Most lines of sight inside the Local Bubble show iron depletions in the range of —0.7 > 
D(Fe) > —1.3, which are typical for warm diffuse clouds. Two lines of sight (36 Oph and 
G191-B2B) show smaller iron depletions and have values of (D/H) gas that are typical for the 
Local Bubble (14.1 ± 5.8 ppm and 16.6 ± 4.1 ppm, respectively). There is no measurement 
of D(Fe) for the line of sight to (3 Cet, but magnesium is undepleted for this line of sight 
(Piskunov et al. 1997). Also, silicon, another refractory metal, is not depleted on this line 
of sight. From Jenkins (2004), one sees that D(Fe) « D(Si) - 0.8, so we might estimate 
D(Fe) w —0.8 for the line of sight to (3 Cet, which also has the highest value of (D/H) gas = 



-19- 



21.9tg;g ppm in the Local Bubble. This connection of high (D/H) gas with minimal metal 
depletion is as predicted by the deuterium-depletion model. 

In Figure 3, we plot (D/H) gas vs. D(Fe) for all 38 lines of sight in Table 3 with mea- 
sured D(Fe) values (excluding the upper and lower limits). Figure 3 shows a clear trend of 
decreasing (D/H) gas with increasing iron depletion, especially for very large Fe depletions, 
D(Fe) < —1.5. The solid line in the plot is the least-squares-weighted linear fit to the data 
as described above. We fit the data with an equation of the form (D/H) gas = (a±l<r) + 
(h±la) D(Fe). For this fit, the parameters are a= 30.00 ± 1.60 and b= 13.04 ± 1.15. The 
Spearman rank correlation test 5 of the unweighted data points shows a highly significant 
correlation parameter and rejects the null hypothesis with 2.9a significance. These param- 
eters are listed in Table 5. While there are five data points that lie above the linear fit, 
the statistical correlation of large Fe depletions with small values of (D/H) gas supports the 
deuterium-depletion hypothesis 6 . In § 6.3 we discuss possible causes of scatter about the 
trend line, which in particular make 7 2 Vel no longer appear anomalous. 

Figure 4 is similar to Figure 3 but with the depletion of silicon D(Si) for the x axis. 
The values of D(Si), listed in column (9) of Table 3, were obtained using iV(Si II) from the 
references and the tabulated values of iV(H I). More significant revisions in iV(Si II) are 
needed than for Fe II as the f-values used for the weak Si II 1020.7 and 1808.0 A lines have 
changed greatly since the Copernicus and IUE eras. We use the Morton (2003) f-values for 
these lines, which are based on laboratory measurements. The changes in iV(Si II) range 
from +0.43 dex for BD+39°3226 to -0.02 dex for /3 CMa with most changes significantly 
positive. We compute the silicon depletion assuming the solar abundance of silicon as log 
(Si/H) = -4.49 ±0.04 (Asplund et al. 2005). We find that D(Si), like D(Fe), shows a clear 
trend of increasing depletion with lower deuterium abundance. Our weighted least-squares 
fit to the 20 data points (excluding upper limits) is of the form (D/H) gas = (a±lcr) + (h±la) 
D(Si). The Spearman rank correlation test of the unweighted data points rejects the null 
hypothesis with 1.4<7 significance. The parameters for this fit are given in Table 5. 

In their recent paper, Prochaska, Tripp, & Howk (2005) find a similar correlation of 
(D/H) gas with the depletion of titanium, which also has a high condensation temperature 



5 In §6.3.5 wc redo this analysis using only the data obtained with the STIS, GHRS, IMAPS, and FUSE 
instruments. The effect is to remove the data points with the largest errors, leaving data points that mostly 
have similar errors and thus similar weights. The result for 26 data points is a similar linear fit that rejects 
the null hypothesis of no correlation with 2.3a significance. 

6 There is no theoretical reason for assuming that the best-fit relation should be linear, or that the D 
depletion and Fe depletion should be tightly correlated with one another. 
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for the lines of sight to seven of the stars listed in Table 2 located beyond the Local Bubble. 
Although Ti II has a similar ionization energy to H I and D I, the velocity profiles of 
the Ti II lines do not track very well the profiles of O I 1335 A for the line of sight to 
HD 195965, perhaps due to different depletion levels as suggested by Prochaska, Tripp, & 
Howk (2005). Although the correlations of (D/H) gas with D(Fe) and D(Si) are significant, the 
scatter that we have seen in Figure 3 and Figure 4 suggests that different physical processes 
are responsible for the depletions of D compared to the metals and most likely among the 
different metals. This could be the cause of the differences between the O I and Ti II line 
profiles for HD 195965 and would suggest significant scatter in plots of (D/H) gas vs D(Ti). 

Wakker & Mathis (2000), Jenkins et al. (1986), and others have demonstrated that the 
depletion of iron and other refractory elements becomes larger with increasing A r (H I), or 
increasing mean hydrogen number density, (%), along the line of sight. This relation is 
generally explained in terms of the increasing contribution of quiescent cold clouds, with 
their large values of iV(H I) and greater metal depletions, to the total iV(H I) when a line of 
sight passes through one or more such cold clouds. See also Oliveira et al. (2006) for a study 
of the correlation of (D/H) gas with (%) along the line of sight. The correlations of (D/H) gas 
with D(Fe) shown in Figure 3 and with D(Si) shown in Figure 4 are consistent with D, Fe, 
and Si being depleted onto grains in cold clouds. 

6.2. Second Test: H 2 Rotational Excitation 

The second empirical test of the deuterium depletion hypothesis is whether or not there is 
a correlation of (D/H) gas with H 2 rotational temperature. Table 6 summarizes the molecular 
hydrogen data for the fourteen lines of sight with (D/H) gas measurements. When both iV(H 2 ) 
and iV(HD) are measured for a given line of sight, we include these molecules in the value 
of (D/H) gas . With some exceptions, for column densities exceeding logiV(H) ps 20.7 (in our 
Galaxy), many H 2 spectral lines become optically thick, and the resulting shielding of the 
FUV radiation field decreases the photodissociation rate leading to a high molecular hydrogen 
fraction, /(H 2 ) (Spitzer & Jenkins 1975; Savage et al. 1977). Under these conditions cloud 
densities can be sufficiently high that the H 2 rotational excitation temperature Tqi, derived 
from the J=l and J=0 levels of H 2 , measures the gas kinetic temperature. In their analysis 
of FUSE spectra of 129 sightlines to O and B stars in the Galactic plane, Shull et al. (2005) 
find that T 01 = 86 ± 20 K. Gillmon et al. (2006) find that T 01 = 124 ± 8 K characterizes lines 
of sight through the halo. For the Galactic plane data set, log [AT (H I) + 2jV(H 2 )] > 20.9, 
while our data set is mostly for smaller column densities. The three lines of sight in Table 6 
with the highest column densities overlap the Shull et al. (2005) data set, and their values 
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of Tqi are consistent with the Shull et al. (2005) range as expected. 

Eleven of the fourteen lines of sight listed in Table 6 have column densities log[A r (H I) + 
2iV(H 2 )] much smaller than 20.9 and, for the most part, values of T 01 much larger than 
86 K. Since this is a new region of parameter space to explore, we decided to test for a 
correlation between T m and (D/H) gas . Given that there is essentially no H 2 located inside 
the Local Bubble, we plot (D/H) gas _ LB vs. T 01 . The data shown in Figure 5 do indeed 
show a correlation with one discrepant point, S Ori, which is located behind a cloud that has 
been excited by stellar winds or explosive events (Jenkins et al. 2000b). Our weighted least- 
squares fit to the 16 data points is of the form (D/H) gas _LB = (a+lo - ) + (bilcr) log (Toi). 
The Spearman rank correlation test of the unweighted data rejects the null hypothesis with 
2.2cr significance, indicating a credible correlation. The fit parameters are listed in Table 5. 
Although the correlation is only at a 2.2a confidence level, we should ask why there should 
be any correlation at all. Since the lines of sight are generally inhomogeneous, most of the 
molecular hydrogen along a given line of sight should be located in cool or cold clouds and 
thus not be cospatial with most of the atomic hydrogen and deuterium. Thus the properties 
of H 2 , such as T m , could be independent of the depletions of D, Fe, and Si averaged over the 
line of sight. Nevertheless, we have found a correlation of T m with (D/H) gas and therefore 
indirectly with the depletions of Fe and Si. Furthermore, for iV(H) < 20.7, T i should 
not be a good measure of the gas temperature but is rather determined by the balance of 
excitations and de-excitations of the lower levels of H 2 with collisions likely unimportant. We 
note that in Table 6 all of the lines of sight with high values of T m have very small fractional 
abundances of H 2 , /(H 2 ). For these lines of sight there is likely very little cold gas and thus 
less depletion of deuterium and metals. We do not know whether this is a good explanation 
for the correlation of (D/H) gas with Tqi shown in Figure 5. Observations of more lines of 
sight with low values of iV(H) but containing H 2 and further theoretical work are needed to 
verify the correlation and to explore other possible explanations. 



6.3. Possible Causes for Scatter and Discrepant Lines of Sight 

Figure 3 shows that five lines of sight (7 2 Vel, Lan 23, WD 1034+001, PG 0038+199, 
and TD1 32709) have high values of (D/H) gas as compared to the general trend of the other 
lines of sight. In Figure 5, the line of sight to 5 Ori has a low value of (D/H) gas _LB compared 
to the trend line. What could explain these discrepant data points and, for that matter, the 
scatter of data about the least-squares trend lines? We list here several considerations that 
may play important roles in explaining the scatter, but this list may not be complete: 
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6.3.1. Inhomogeneous Lines of Sight 

Until now, we have assumed that the properties of interstellar gas along a given line 
of sight have constant values inside the Local Bubble and then a different set of constant 
values beyond the Local Bubble. This very simple approximation is invalid for most lines 
of sight that have been studied with sufficient spectral resolution to identify the different 
velocity components and their individual properties. If the properties of interstellar gas in 
the different components are significantly different, for example in temperature, ionization 
(H II vs. H I regions), and density, then different observables may be formed preferentially in 
different locations and thus not correlate well. One example is the correlation of (D/H) gas _LB 
with the temperature of H 2 , as H 2 could be present in only one of the components along 
the lines of sight. For lines of sight with low H2 column densities, H2 could be found in 
relatively warm gas associated with H I and neutral metals. For example, high-resolution 
spectra of 5 Ori obtained with the IMAPS instrument show H 2 line profiles with similar 
velocity structure as N I, O I, and D I (Jenkins et al. 1999, 2000b). This could explain why 
T i is large for the the 5 Ori line of sight. 

One way of estimating the depletion of an element in a particular velocity component 
without information on N(H. I) for that component is to reference the gas-phase column 
density of the element to another element like sulphur, which is not significantly depleted 
in H I regions. Howk et al. (1999) have used this approach to derive D(Fe) and D(Si) for 
Component 1 in the line of sight to fi Col using GHRS echelle spectra. This low-velocity 
component, also called Component A by Shull & York (1977), contains most of the neutral 
gas and H 2 along the line of sight. Tables 2, 3, and 6 include the (Howk et al. 1999) data 
for Component 1. 



6.3.2. Line Saturation 

For complex lines of sight with large iV(H I), it is important to measure iV(D I) from the 
higher lines in the Lyman series to avoid line saturation that could lead to underestimates in 
iV(D I) and therefore in (D/H) gas . The 20 km s _1 spectral resolution of FUSE is inadequate 
to resolve the D lines or to separate velocity components in the line. In particular, narrow 
components from cold gas could be saturated but not recognizable at the resolution of FUSE. 
We have argued in § 5 that the low values of (D/H) gas for lines of sight beyond the Local 
Bubble could be due to the dominant role of cool diffuse clouds in which deuterium is 
largely depleted. Such clouds would have narrow absorption features observable in metal 
lines. Indeed, Redfield & Linsky (2002) list a number of lines of sight with Ca II absorption 
features with b < 1.5 kms -1 . The chances of this occuring grow with increasing N(R I), 
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but measurements of iV(D I) in high members of the Lyman series, which require high 
S/N data, decrease the probability that line saturation will occur. One can test for D line 
saturation by looking for an increase in N(D I) toward higher lines in the Lyman series or 
from a curve of growth analysis (e.g., Friedman et al. 2002). Also, high resolution spectra of 
interstellar metal lines can identify velocity features that could have saturated D absorption. 
The blending of adjacent Lyman lines with each other and with the many H2 lines will 
eventually set a limit to the largest value of 7V(H I) for which one can measure unsaturated 
D Lyman lines, but this should not occur until logiV(H I) is considerably larger than 21.0 
(Jenkins 1996). 



6.3.3. Ionization Corrections 

Jenkins et al. (1986) and Jenkins (2004) have argued that anomalously small metal 
depletions compared to typical lines of sight with the same N(H I) could be due to H II 
regions where hydrogen is partially ionized but metals with second ionization potentials 
greater than 13.58 eV, such as iron, silicon, and phosphorus (Lehner et al. 2004), remain 
mostly singly ionized. Positive values of D(Si) are consistent with Si + being the dominant 
ionization state and H being partially ionized. A good example of this is the line of sight to 
/3 CMa. If we assume that essentially all of the hydrogen is H I and all of the silicon and iron 
are singly ionized for this line of sight, then D(Si)= +0.41 ±0.16, and D(Fe) = —0.30 ±0.16 
(Gry, York, & Vidal-Madjar 1985; Dupin & Gry 1998). However, along this line of sight 
hydrogen is mostly ionized (Gry, York, & Vidal-Madjar 1985). In their analysis of this line 
of sight, Dupin & Gry (1998) identify four velocity components with most of the H I column 
density in Component C. We therefore list in Table 3 the depletions for Component C that 
include all ionization stages for hydrogen, silicon, and iron, D(Si) = ±0.041°;^ and D(Fe) 
= — 0.61^q ' 14 (Dupin & Gry 1998). The line of sight to (3 Cet also has a positive value for the 
silicon depletion, D(Si) = ±0.30 ± 0.41. Redfield & Linsky (2004b) showed that one of the 
two velocity components toward /3 Cet has a very high temperature, T = 12, 400 ± 2, 800 K, 
likely indicating partial ionization of hydrogen. Thus ionization corrections of presently 
unknown size should be applied to D(Si) for this line of sight, which would move its data 
point to the right in Figure 4. The high temperature of the gas toward (3 Cet is likely due to 
a recent shock, which is consistent with the high value of (D/H) gas = 21.9+lg ppm for this 
line of sight. 

The 7 2 Vel line of sight has seven velocity components - three are H II regions, likely 
produced by ionization from the Wolf-Rayet star in the binary system, and four are H I 
regions (Fitzpatrick & Spitzer 1994). The H I regions contribute 89% of iV(Fe II), but the 
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fractional ionization of hydrogen and iV(H I) in each component is not known. The iron 
depletion listed in Table 3 was therefore computed from the line-of-sight integrated values 
of 7V(H I) and N(Fe II), rather than from the column densities in the H I components. 
Inclusion of the unknown iV(H II) would shift D(Fe) to the right in Figure 3. Since sulfur is 
not usually depleted, we estimate that 72% of the total hydrogen column density is located 
in the H II regions by summing the column densities of S I, S II, and S III in each velocity 
component. Although the inclusion of this amount of ionized hydrogen would shift D(Fe) for 
7 2 Vel by 0.55 dex to the right in Figure 3, an unknown but considerable amount of Fe could 
be doubly ionized, shifting D(Fe) by an unknown amount to the left. Realistic calculations 
of the ionization of H and Fe are needed to address this problem. 

Ionization corrections can be important for lines of sight with small values of iV(H I) 
because ionizing radiation can then penetrate the gas and preferentially ionize H compared to 
Fe + and Si + . For example, Frisch & Slavin (2003) find in their ionization equilibrium models 
for the nearby ISM that hydrogen is about 30% ionized while 96.5% of Fe is Fe II and 99.6% 
of Si is Si II. To avoid lines of sight where such ionization corrections are likely important, 
we reconsider the correlation of (D/H) gas with D(Fe) and D(Si) but now remove the data 
for lines of sight with logiV(H I) < 19.0. This removes all but two of the Local Bubble lines 
of sight. In Figure 6, we plot (D/H) gas vs. D(Fe) for the 24 data points that meet this 
criterion. The slope is now steeper than for the full data set (see Table 5). The Spearman 
rank correlation test rejects the null hypothesis of no correlation with 2.6cr significance. 

Figure 7 is similar to Figure 4, except that we remove the data points with log iV(H I) < 
19.0 and plot (D/H) gas vs. D(Si) for the 11 data points that meet this criterion. The slope 
is nearly the same as for the full data set (see Table 5). The Spearman rank correlation test 
rejects the null hypothesis of no correlation with 2.1<r significance. The data are in excellent 
agreement with the linear fit, except for the one uncertain data point, /i Col. 

The removal of lines of sight that are most likely affected by ionization corrections leads 
to steeper correlations of (D/H) gas with D(Fe) and D(Si), fewer discrepant points, and tighter 
correlations. The remaining scatter could be due to measurement errors (especially for some 
of the Copernicus data), inhomogeneous lines of sight, and different grain compositions. The 
ionization of H and D should be the same as they share the same bound-free continuum. 

6.3.4- Different Grain Compositions 

One should not expect a simple one-to-one relation between (D/H) gas and D(Fe) or 
D(Si) along a given line of sight, because deuterium and different metals may be depleted 
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onto different types of grains or be located preferentially in different layers of the same grain. 
Many authors (e.g., Spitzer & Fitzpatrick 1993) have argued that grains typically have hardy 
cores but more easily destroyed mantles. For example, iron may be concentrated in the grain 
cores and some deuterium concentrated in PADs either free-flying or incorporated into larger 
grains. Jones et al. (1994) estimated the destruction timescale for grains by supernova-driven 
shocks to be 2-3 xlO 8 yr but the timescale for the erosion of grain mantles by weaker shocks 
is much shorter. We note that the five lines of sight that lie above the trend line in Figure 3 
(7 2 Vel, Lan 23, WD 1034+001, PG 0038+199, Tl 32709) also have the highest values of 
(D/H) gas (see Figure 1). (Feige 110, the fourth line of sight with a high value of (D/H) gas , does 
not have D(Fe) or D(Si) measurements.) Clearly, there is no simple one-to-one relationship 
between (D/H) gas and metal depletion, although there are the general trends indicated by 
the least-squares trend lines. In Figure 6, where we have considered only lines of sight with 
logiV(H I) > 19.0 to minimize ionization corrections, all of the five high lines of sight are 
still present. 

In Figure 3 the dashed line is drawn parallel to the least-squares fit to all of the lines of 
sight (solid line) but arbitrarily displaced upwards by 8.5 ppm. The dashed line is a good fit 
to the five high data points. The high values of (D/H) gas for these five lines of sight could be 
explained by either a smaller than usual percentage of the grains being composed of carbon 
and thus few sites available for deuterium to deplete on or the evaporation of deuterium 
that was in grain mantles by weak shocks. If the slope of the dashed line is real and not an 
artifact of the few data points, then there is support for a smaller than usual percentage of 
carbon grains in these lines of sight. Since the condensation temperature for carbon is much 
smaller than that for iron (typically in the form of oxides) (Sofia, Cardelli, & Savage 1994; 
Savage & Sembach 1996), a low percentage of carbon grains compared to iron grains is likely 
for some lines of sight. 



6.3.5. Difficulty of Measuring N(H I) 

The most difficult column density to measure accurately is usually hydrogen because 
the Lyman lines are very optically thick, they must be measured against an uncertain stellar 
background (Lyman absorption or emission lines), and they are very broad, which requires 
interpolation over a large wavelength range. The hydrogen column densities cited in Table 2 
are based on careful analyses that consider these and other effects. The most accurate 
values of iV(H I) are generally obtained from analyses of STIS, GHRS, and FUSE spectra. 
We note that iV(H I) for the Feige 110 line of sight was obtained from the analysis of one 
IV E spectrum. Since the Lyman-a line lies near the end of the IV E spectrum where the 
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sensitivity is low and the echelle spectral orders are close together, it is difficult measure the 
background and continuum. While there is no evidence that the value of iV(H I) listed in 
Table 2 is more uncertain than the cited error bars, the Feige 110 line of sight should be 
reobserved with a better spectrograph when feasible. The Lan 23 line of sight should also be 
reobserved as iV(H I) was obtained from a noisy low dispersion EUVE spectrum rather than 
from Lyman line spectra. To test whether the less reliable values of iV(H I) are biasing our 
conclusions, we plot in Figure 8 only the data points for which iV(H I) was obtained from 
STIS, GHRS, or FUSE. The fit to the data and the parameters characterizing the fit listed 
in Table 5 are similar to the results obtained when we included all of the data in Figure 3. 



6. 3. 6. Removing a Possible Correlation from the Plots 

In Figures 3-4 and 6-8, the measured 7V(H I) values enter the quantities plotted in both 
the x and y-axes, raising the possibility of inaccurate or false correlations. To eliminate this 
possibility, we have plotted in Figure 9 the log[iV(D I)/JV(Fe II)] ratio vs. logiV(H I) for all 
lines of sight. We include in the figure the weighted least-squares fit to all of the data points 
(solid line) and the fit to only the lines of sight which extend beyond the Local Bubble (dashed 
line). The essentially zero slope in the fit to the data beyond the Local Bubble indicates 
that the positive slope of the fit to the entire data set is determined by the low values of 
log[iV(D I)/JV(Fe II)] for lines of sight inside the Local Bubble. The zero slope beyond the 
Local Bubble indicates that on a statistical basis the deuterium and iron depletions are 
proportional. The scatter of the data about the two regression lines suggests different grain 
compositions along the various lines of sight or different processes that deplete deuterium 
and iron from the gas phase but that depletion is important for both species. While the 
scatter indicates that the depletions of D and Fe are not proportional on individual lines of 
sight, the zero slope fit to the data beyond the Local Bubble indicates that on average the 
rates of depletion of D and Fe appear to be proportional. 

We noted in §6.3.3 that in regions of low shielding (logiV(H I) < 19.0) from Lyman 
continuum radiation primarily from hot stars, that hydrogen can be partially ionized, whereas 
metals with second ionization potentials > 13.58 eV will be mostly singly ionized. The 
second ionization potential of iron is 16.18 eV, whereas the ionization potential of neutral 
deuterium is the same as that of hydrogen. Thus partial ionization of deuterium could explain 
the decrease in N(D I)/iV(Fe II) at logiV(H I) < 19.0. It may also explain the very low 
ratios for the short lines of sight to 36 Oph and G191-B2B. Redfield & Linsky (2000) found 
that hydrogen is half ionized in the Local Interstellar Cloud, which has a maximum column 
density logiV(H I) = 18.3. This result could simply explain the factor of two decrease in 
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N(D I)/iV(Fe II) between lines of sight with large and small values of iV(H I) in Figure 9. 

7. WHAT IS THE TOTAL DEUTERIUM ABUNDANCE IN THE LOCAL 

GALACTIC DISK? 

We propose that the most plausible explanation for the wide range in (D/H) gas mea- 
surements within about 1 kpc of the Sun is that different amounts of deuterium depletion 
occur in different lines of sight. We have come to this conclusion because: (i) theoretical 
arguments suggest that deuterium atoms can replace enough hydrogen atoms in carbona- 
ceous grain materials to explain low values of (D/H) gas _ LB in the undisturbed ISM, (ii) high 
(D/H)dust ratios have been measured in interplanetary dust particles that likely came from 
the ISM, providing an important proof of concept for the deuterium depletion hypothesis, 
and (iii) the correlation of large iron and silicon depletions in the gas phase, D(Fe) and 
D(Si), with low values of (D/H) gas strongly support the hypothesis that dust grains are an 
important reservoir for deuterium and that deuterium returns to the gas phase when the 
grains are disturbed. 

The deuterium-depletion model predicts that the most likely value for the ratio of deu- 
terium in all forms (gas plus dust) to hydrogen, (D/H) to tai, for the local region of the Galactic 
disk, (D/H) LDtot , should be equal to or slightly above the highest measured (D/H) gas ratios. 
If some deuterium remains on grains for these lines of sight, then the total amount of deu- 
terium will slightly exceed the gas-phase value. In Figure 1 there are five lines of sight (7 2 Vel, 
Lan 23, WD 1034+001, Feige 110, and LSE 44) extending beyond the Local Bubble that 
have the highest values of (D/H) gas 22 ppm. These are the best candidates for inferring 
(D/H)LD tot in the local Galactic disk. The weighted mean and standard deviation of the 
mean for these five data points is ((D/H) gas ) = 21.7 ± 1.7 ppm. However, the Local Bubble 
foreground has a known mean value of (D/H) gas and extent, logiV(H I) ps 19.2, indicating 
that deuterium is depleted inside the Local Bubble. We therefore subtract the Local Bubble 
foreground from the (D/H) gas values to obtain the values of (D/H) gas _ LB listed in Table 3. 
The weighted mean values for these five lines of sight is ((D/H) gas _ L B) = 23.7 ± 2.4 ppm. 
The line of sight to 7 2 Vel lies closest to the Local Bubble (log N(H I) = 19.710 ±0.026) and 
is thus most subject to the systematic errors (especially the uncertain extent of the Local 
Bubble in iV(H I)) associated with subtracting the Local Bubble contribution. However, 
this line of sight has the highest weight (58%) in computing this mean. To minimize this 
potential source of systematic error, we have recomputed the mean without weighting by 
the inverse errors, obtaining ((D/H) gas _ L B) = 23.1 ± 2.4 ppm. Given that we have selected 
the data points to be used in computing the mean from a larger sample, the errors are not 
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Gaussian, but we believe that our method gives the most representative value (D/H) gas _LB 
available at this time. Since even for these five lines of sight some deuterium could be in the 
grains, (D/H) LDtot > ((D/H) gas _ LB ) = 23.1 ± 2.4 ppm. 

Even the higher estimate of the primordial D/H ratio (see Section 1), (D/H) prim = 
27.5^9 ppm, implies a small deuterium astration (or depletion) factor for the local region of 
the local Galactic disk, f d =(D/H) prim /(D/H) LDtot < (27.5+^)/(23.1 ± 2.4) < 1.19+°;^. For 
the lower estimate of (D/H) prim = 26.01J;? ppm, f d < (26.0l^)/(23.1 ± 2.4) < 1.12 ± 0.14 
These small values of fd provide an important constraint on models of Galactic chemical 
evolution, which must also explain the evolution of many chemical species. The models 
of Chiappini, Renda, & Matteucci (2002) (cf. Romano et al. 2003), for example, predict 
fd ~ 1.5 for a wide range of assumed (D/H) pr j m and different rates of extra mixing in stellar 
interiors. These models include two epochs of infall of D-rich and metal-poor gas 7 from the 
halo or IGM, and the models predict sensible present epoch values for 3 He/H, 4 He/H, and 
7 Li/H. By altering the assumptions regarding the rates of infall and star formation, more 
recent work by Romano et al. (2006) is able to reproduce current data on stellar abundance 
patterns in the solar neighborhood (R = 8 kpc) with astration factors 1.3 < fd < 1-8; the 
lower end of this range overlaps our determinations fd < 1.19^^1 or f d < 1.12 ± 0.14. 

There are also estimates of the D/H ratio in the protosolar cloud, (D/H) psc , when the 
Galaxy was about two-thirds of its present age. Geiss & Gloeckler (1998) have inferred 
(D/H) psc = 21 ±5 ppm based on measurements of 3 He/ 4 He in the solar wind, measurements 
of the same quantity in Jupiter's atmsphere, and the assumption that 3 He in the solar 
atmosphere has not undergone nuclear reactions. Analysis of Infrared Space Observatory 
(ISO) spectra of H 2 , HD, CH 4 , and CH 3 D in Jupiter's atmosphere led Lellouch et al. (2001) 
to infer (D/H) psc = 21 ± 4 ppm, although there are many assumptions that go into this 
analysis. One cannot directly compare (D/H) psc with present day estimates of (D/H) L Dtot, 
because the Sun has moved a considerable distance from its place of birth about 4.6 Gyr 
ago. However, the similar values of (D/H) psc and (D/H)LDtot > 23.1 ± 2.4 ppm suggest that 
most of the decrease in the total deuterium abundance in the ISM may have occurred as a 
result of nuclear reactions in stars and subsequent transfer of deuterium-depleted gas into 
the ISM during the first two-thirds of the age of the Galaxy. 



7 Thc infalling gas has a metal abundance 20% solar (Romano et al. 2003), which is consistent with the 
results for the high velocity cloud Complex C (Collins et al. 2003). 
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8. COULD INFALL EXPLAIN THE HIGH (D/H) gas RATIOS? 

It is fair to ask whether depletion is the only operative mechanism for changing the 
abundance of deuterium. Oliveira et al. (2005) have used the constancy of the (0/H) gas 
ratios over a range of environments in the local part of the Galaxy to argue that significant 
changes in the D /H ratios because of infall of lower metallicity gas are unlikely While it is 
tempting to suggest significant and highly localized infall for the high D/H sightlines, such 
a large infall would dramatically lower the O I abundances along these sightlines. A large 
increase in the value of (D/H) gas from 15 ppm to 22 ppm would require a highly localized 
infall of gas with (D/H) tota j = 24 ppm within a mixing time of about 350 Myr (de Avillez & 
Mac Low 2002). The fraction due to infall would equal more than three quarters of the mass 
of the gas along a given sightline. Such a large infall of D-rich, but metal-poor, gas would 
also lower the measured (0/H) gas ratios by more than a factor of three; such low values are 
not observed. The constancy of (0/H) gas compared to variations in (D/H) gas thus argues 
against localized infall (Oliveira et al. 2005). Although the infall mechanism is unlikely to 
explain the wide range of observed (D/H) gas values, we do not rule out the possibility of 
another, as yet unknown, mechanism for increasing the D abundance. 

Finally, we include a note of caution. Hebrard et al. (2005) have pointed out that the 
(D/0) gas data show variations similar to those observed for (D/H) gas , except that there were 
no high values at that time corresponding to those measured for (D/H) gas . The relative 
uniformity of O I/H I allows O I to be used as a proxy for H I, with the advantage that 
for a given line of sight, N(0 I) is much closer to iV(D I) than iV(H I). Therefore, the D I 
and the O I transitions in the FUSE range have similar opacities, lowering the chances of 
systematic errors due to large differences in optical depth (Hebrard & Moos 2003). A gas- 
phase (D/H) gas ratio of 23.1 ± 2.4 ppm and the (0/H) gas ratio of 343 ± 15 ppm measured by 
Meyer et al. (1998) implies (D/0) gas = 7.02 x 10 -2 , nearly twice the mean measured value of 
(3.84 ±0.16) x 10~ 2 (Hebrard & Moos 2003). Although such high values had not previously 
been published, Oliveira et al. (2006) report 3 new sightlines (WD 1034+001, BD+39 3226, 
and TD1 32709) with high values of (D/0) gas , consistent with (D/0) gas « 7.02 x 10~ 2 . 

9. CONCLUSIONS 

The main goal of the FUSE mission has been to obtain accurate measurements of 
(D/H) gas for many sightlines in the Milky Way Galaxy and beyond in order to measure 
(D/H) prim and to obtain constraints on Galactic chemical evolution. FUSEhas now obtained 
the data needed for accurate measurements of (D/H) gas much further out in the Galaxy 
and with higher precision than previous instruments. In our summary of the published 
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measurements of (D/H) gas obtained from FUSE, HST, Copernicus, and IMAPS, we find, 
following Wood et al. (2004), that the (D/H) gas measurements appear to fall into three 
distinct groups depending on the neutral hydrogen column density iV(H I) to the target 
star. For lines of sight with logiV(H I) < 19.2, that is within the Local Bubble, the mean 
value of (D/H) gas is 15.6 ± 0.4 ppm, where the uncertainty is the standard deviation of the 
mean. At large column densities, logiV(H I) > 20.7, we find that the five lines of sight have 
much lower values of (D/H) gas with a mean value of 8.6 ±0.8 ppm. In the intermediate regime 
there is a factor of 4-5 range in the high precision (D/H) sas measurements with values as 
low as 5. Olg^ ppm (6 Car) and as high as 21.8 ± 2.1 ppm (7 2 Vel). Since this large range 
has been measured from spectra obtained by several instruments, it is not an instrumental 
artifact, but rather an observational result that requires an explantion. 

Both large changes in the local astration rates and infall of D-rich and metal-poor gas 
from the halo or IGM cannot explain the large range in (D/H) gas without predicting large 
variations in the local ISM metal abundances, in particular (0/H) gas , which are not ob- 
served. Instead, we describe a time-dependent deuterium-depletion model in which much 
of the deuterium in the ISM can reside on grains under the appropriate conditions. The 
theoretical basis for this model is the lower zero-point energy of CD bonds compared to 
C-H bonds. When the ISM is undisturbed for a long time, grains are very cold, and deu- 
terium can replace hydrogen on carbonaceous grains and perhaps other types of grains. We 
estimate that D/H in the grains can approach ~0.2, thereby reducing (D/H) gas to 6 ppm or 
lower. When a region of the ISM is shocked or the gas comes close to a hot star, the grains 
are partially or fully destroyed, sending deuterium atoms back into the gas phase. In this 
model, the value of (D/H) gas for a given line of sight depends on the environment and past 
history of the grains. 

Evidence in support of this model includes the measurement of high D/H ratios in 
interplanetary carbonaceous dust grains, which are likely interstellar in origin. Since the 
conditions required from deuterium depletion are similar to those for depletion of metals like 
Fe, Si, and Ti in the ISM, the model predicts a correlation of depletions of these metals with 
low values of (D/H) gas . We present data that strongly support the correlation for Fe and 
Si, and Prochaska, Tripp, & Howk (2005) present the data in support of correlation with 
Ti. The model also predicts that (D/H) gas should be correlated with the recent thermal and 
ionization history of the gas, whereas the variable astration hypothesis would not be expected 
to show such a connection; we find that (D/H) gas has a significant positive correlation with 
the rotational temperature of the molecular hydrogen on the line-of-sight, further supporting 
the hypothesis that (D/H) gas is affected by interstellar processes acting on relatively short 
time scales. We note that in the most recent discussion of Galactic chemical evolution 
models, Romano et al. (2006) conclude that "depletion of deuterium on to dust grains is the 
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most likely physical mechanism proposed so far to explain the observed dispersion in the 
local data." 

If this model is realistic, then the best estimate of the total deuterium abundance, 
(D/H) tota i would be the highest measured values of (D/H) gas because these would be for 
lines of sight with minimum depletion of deuterium. We note that there are five lines of 
sight in the intermediate iV(H I) regime with high values of (D/H) gas . The weighted mean 
value for these four lines of sight is 21.7 ±1.7 ppm (standard deviation of the mean). Since 
there could be some deuterium depletion even in these lines of sight, we conclude that the 
best estimate of the total D/H in the local disk region of the Galaxy, after subtracting the 
Local Bubble foreground column densities, is (D/H) LDto t > 23.1 ± 2.4 ppm 8 . This new 
estimate is > 58% higher than the value of 15 ppm recently used in some Galactic chemical 
evolution models (e.g., Romano et al. 2003), which is (D/H) gas for the Local Bubble rather 
than (D/H) LD tot- 

Two methods for determining (D/H) prim based on data from very different times in the 
early universe are now in agreement. Measurements of (D/H) gas in five quasar absorption 
line systems have a mean value of 27.8^3;g ppm (Kirkman et al. 2003). Analysis of data 
from the WMAP and other cosmic microwave experiments yield Vt b h 2 = 0.0224. Depending 
on the adopted nuclear reaction rates, (D/H) pr i m = 27.5^'g ppm (Cyburt, Fields, & Olive 
2003) or 26.01J;? PP m ( Coc et aL 2004 )- 

Galactic chemical evolution models attempt to explain the decrease in D/H from the 
protogalaxy, for which D/H is presumed to be the primordial value, to the present epoch 
where D/H = (D/H) LDtot- Nuclear reactions in stars and subsequent return of D-poor and 
metal-rich gas by SN explosions and stellar winds will reduce deuterium in the ISM by the 
factor f d = (D/H) prim /(D/H) LDtot . Models of Chiappini, Renda, & Matteucci (2002), for 
example, predict that fd ~ 1.5 for a wide range of assumptions concerning the value of 
(D/H) pr j m and amount of mixing in stellar interiors. This result can be compared with our 
results, f d < (27.5lf.g)/(23.1 ± 2.4) < 1.19lg;^ or f d < (26.0l^)/(23.1 ± 2.4) < 1.12 ± 
0.14. These two results differ by > 1.5a and > 1.8a, respectively, from the predictions of 
the Chiappini, Renda, & Matteucci (2002) models. If, on the other hand, the low value 



8 Very recently Rogers et al. (2005) detected the hyperfine ground-state transition of atomic deuterium 
at 327 MHz (92 cm) in the Galactic anticenter direction. They find that (D/H) gas = 23 ± 4 ppm (la) or 
23+U PP m (3c), and estimate that the column density- weighted mean distance of the emitting gas is about 
2 kpc. Since the beamwidth is 14°, the beam averages an unknown amount of cloud and diffuse gas in and 
above the Galactic plane. While it is difficult to compare this new result with our line of sight measurements 
over shorter paths, the value of (D/H) gas obtained from the hyperfine transitions of D and H is consistent 
with our proposed value of (D/H) LDtot- 
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of (D/H) gas = 8.6 ± 0.8 ppm for the five most distant lines of sight is representative of 
(D/H) LDtot as Hebrard et al. (2005) proposes, then f d would be (27.5j£|)/(8.6 ± 0.8) = 
3. 21^37 or (26.0lJ;?)/(8.6±0.8) = 3.0 ±0.35, which are far larger values than the Chiappini, 
Renda, & Matteucci (2002) models predict. Very recent Galactic chemical evolution models 
(Romano et al. 2006) can now accomodate values of fd as small as ~1.3, consistent with our 
determination fd < LlQ+aijj or 1-12 ± 0.14, depending on the adopted value of (D/H) prim . 
However, our empirical results for f d are smaller than f d = 1.39, the lowest value for a 
model discussed by Romano et al. (2006). We suggest that new Galactic chemical evolution 
models should be developed to see whether smaller values of fd are consistent with the other 
empirical abundance and theoretical constraints that the models must fit. 

We also call attention to the need for more high quality data to better constrain 
(D/H) L Dtot and (D/H) prim . In particular, we need better quality observations of the highly 
saturated Lyman lines for which measurements of N(R I) may contain systematic errors 
of an unknown nature. Also, the assumptions underlying the Galactic chemical evolution 
models should be reexamined to see whether fd might be as small as ~1.25 without violating 
other constraints. 
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Table I. PRIMORDIAL D/H: MEASURED AND INFERRED a 



References 


Q B h 2 
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(D/H) prim 


Data Used b 








(ppm) 





Kirkman et al. (2003) 0.0214 ± 0.0020 5.9 ± 0.5 27.8±;j;8 5 quasar absorption line systems 

Spergel et al. (2003) 0.0224 ± 0.0009 6.1±g;l 26.2l^ CMB, 2dFGRS, and La forest data 

Sanchez et al. (2006) 0.0225 ± 0.0010 6.15 ± 0.27 25.9±^ CMB, 2dFGRS, 6 parameter model 

Cyburt, Fields, & Olive (2003) 0.0224 ± 0.0009 6.14 ± 0.25 27.5±H CMB, new nuclear reaction rates c 

Coc et al. (2004) 0.0224 ± 0.0009 6.14 ± 0.25 26.0±};? WMAP, SBBN, new reaction rates d 



a Quoted uncertainties are ±lcx. 

b CMB refers to the cosmic microwave background data consisting primarily of data from the Wilkinson Microwave 
Anisotropy Probe (WMAP) but also finer scale CMB experiments. 2dFGRS refers to the power spectrum of galaxy 
clustering measured from the final two-degree field galaxy redshift survey (Sanchez et al. 2006). SBBN refers to 
standard big band nucleosynthesis. 

c Nuclear reaction rates of Angulo et al. (1999). 

d Nuclear reaction rates cited in Coc et al. (2004). 
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5,27,28 


G191-B2B 


156 


+07 


69 ± 15 


18.18 ± 0.09 


13.40 ± 0.035 


FUSE 


LBg 


13.05 ± 0.02 


12.55+°'li 

— 0.14 


9,27,28 


(3 CMa c 


226 


-14 


153 ± 15 


18.20+°-" 
u w — 0.20 


> 13.40 


Copernicus 


LBg 


13.04 ± 0.04 


13.745 ± 0.04 


16,29 


(j Gem 


191 


+23 


37 + 1 


18.201 ± 0.037 


13.34 ± 0.05 


HST 


LBg 






1 


Capella 


163 


+05 


12.9 ± 0.1 


18.239 ± 0.035 


13.44 ± 0.02 


HST 


LBg 


12.49 ± 0.02 


13.00 ± 0.05 


2,27,28 


f3 Gem 


192 


+23 


10.34 ± 0.09 


18.261 + 0.037 


13.43 + 0.05 


HST 


LBg 


12.42 + 0.03 




1,27 


a Tri 


139 


-31 


19.7 + 0.3 


18.327 + 0.035 


13.45 + 0.05 


HST 


LBg 


12.65 + 0.20 




1,27 


(3 Cct 


111 


-81 


29.4 + 0.7 


18.36 + 0.05 


13.701°°? 


HST 


LBg 




14.17 + 0.41 


23,28 


A And 


110 


-15 


25.8 + 0.5 


18.45 ± 0.075 


13.68 + 0.10 


HST 


LBg 






3 


Feige 24 


166 


-50 


74 + 20 


18.47 + 0.03 


13.58 + 0.14 


HST 


LBg 




13.54 + 0.07 


10 


WD 0621-376 


245 


-21 


78 ± 23 


18.70 + 0.15 


13.85 ± 0.045 


FUSE 


LBg 


12 94+° 14 
rz.yi_ 16 




11 


WD 2211-495 


346 


-53 


53 ± 16 


18.76 + 0.15 


13.94 + 0.05 


FUSE 


LBg 


13.25 + 0.075 


14.00 + 0.08 


7,30 


WD 1634-573 


330 


-07 


37 + 3 


18.85 ± 0.06 


14.05 ± 0.025 


FUSE 


LBg 


13.08 + 0.065 


13.761°;" 


6,30 


a Vir 


316 


+51 


80 + 6 


19.0 + 0.1 


14 2+°' 2 


Copernicus 


LBg 


13 16+ 001 


I4.l6t°;° 8 


13,32 


GD 246 


87 


-45 


79 + 24 


19.11 ± 0.025 


14.29 ± 0.045 


FUSE 


LBg 


13.30 + 0.05 


14 07+0.06 


12 


A Sco d 


352 


-02 


216 ± 42 


19.23 + 0.03 


14 11+0.09 


Copernicus 


INT 


13.04 + 0.04 


13.42 + 0.04 


19 


(3 Cen 


312 


+01 


161 ± 15 


19.63 + 0.10 


14.7 + 0.2 


IUE, Copernicus 


INT 


13 09+0 05 
io.yz_ ()4 


14.59 + 0.01 


13,32 


7 2 Vel 


263 


-08 


258 ± 35 


19.710 + 0.026 


15.05 + 0.03 


IMAPS 


INT 


13.95 + 0.11 


15.03 + 0.04 


20,26 


a Cm 


300 


+00 


98 + 6 


i9.85tg : ;; 


14.95 + 0.05 


Copernicus 


INT 


14.00 + 0.10 




13,37 


BD+28°4211 


82 


-19 


104 ± 18 


19.846 + 0.018 


14.99 ± 0.025 


FUSE 


INT 


14.10 + 0.10 




14,38,39 


Lan 23 


108 


-01 


122 ± 37 


19.89tg;g 


15.23 + 0.065 


FUSE, EUVE 


INT 


14.03 + 0.065 




12,41 


It Col c 


237 


-27 


400t 7 °° 


19.86 + 0.015 


14.71°;? 


HST, Coper. 


INT 


14.13 + 0.02 


15.10 + 0.02 


13,35 


(PuP 


256 


-05 


429 ± 94 


19.963 ± 0.026 


15.11 + 0.06 


IMAPS 


INT 


U l-j+0.17 


15 07+ 07 
ia.u/_ 02 


20,32 


TD1 32709 


233 


+28 


520 ± 90 


20.03 + 0.10 


15.30 + 0.05 


FUSE 


INT 


13.95 + 0.10 




44 


WD1 1034+001 


248 


+48 


1551" 


20.07 + 0.07 


15.40 + 0.07 


FUSE 


INT 


14.10 + 0.10 




44 



Table 2 — Continued 



Target 


1 

(dcg) 


b 

(deg) 


d 

(PC) 


logJV(H I) a 


logiV(D I) a 


Satellite 


Flag b 


log N(Fe II) a 


logJV(Si II) a 


Refs. 


BD+39°3226 


65 


+29 




20.08 + 0.09 


15.15 + 0.05 


FUSE 


INT 


14.15 + 0.07 


14.80 ±0.20 


31,44 


Feige 110 


74 


-59 


179+^f 


20.14+"'* 5 


15.47 + 0.03 


FUSE, IUE 


INT 






17,38,39 


7 Cas 


124 


-02 


188 ± 20 


20 16+ o: ° 8 


io.io_ 05 


Copernicus 


INT 


14.08tg;S 




18,37 


i Ori 


210 


-20 


407 ± 127 


20 - 15 -007 


-i r + 04 

15-30^55 


Copernicus 


INT 


_, „ „ _i_Q 2() 


-i i- 1 02 

15 - 16 Iao3 


22,32,37,38 


5 Ori 


204 


-18 


281 ± 65 


20.19 + 0.03 


15 06+ 06 
lO.UD_ Q4 


IMAPS 


INT 


14.08 ± 0.03 




21,32,40 


8 Car f 


290 


-05 


135 + 9 


20.28 ± 0.10 


14 QS+0' 18 
141 - a °-0 .21 


Copernicus 


INT 


14.19 + 0.03 


14.38±°;£j 


15 


e Ori 


205 


-17 


412 ± 154 


20 45+0-08 


15.25 + 0.05 


Copernicus 


INT 


14.20 + 0.10 


15.02±°;°« 


22,32,37 


PG 0038+199 g 


120 


-43 


9q7 +l64 


20.48 ± 0.04 


15.76 + 0.04 


FUSE 


INT 


14 42+0- 03 

±4 -* z -0.02 




25 


LSE 44 


313 


+13 


554 ± 66 


20.52±°;i° 


15.87 + 0.04 


FUSE, IUE 


INT 




43 


JL 9 


323 


-27 


590 ± 160 


20.78 ± 0.05 


15.78 + 0.06 


FUSE 


LDg 


14.69 ± 0.085 




23 


HD 195965 


86 


+05 


794 ± 200 


20.95 ± 0.025 


15.88 + 0.07 


FUSE 


LDg 


14.81 ±0.01 




24,38 


LSS 1274 


277 


-05 


580 ± 100 


20.98 ± 0.04 


15.86 + 0.09 


FUSE 


LDg 


14.81 ±0.075 




23,39 


HD 191877 


62 


-06 


2200 ± 550 


21.05 ± 0.05 


15 94+o.n 
ro.»i_ 06 


FUSE 


LDg 


14.95 ± 0.02 




24,38 


HD 90087 g 


285 


-02 


2740 ± 800 


21.22 ± 0.05 


16.16 + 0.06 


FUSE 


LDg 


15.22 ±0.02 




41 



a Quoted uncertainties assumed to be la (see text). iV(H I) is the most recent published value. N(D I) is not listed when not in original paper. 
iV(Fc II) and ./V(Si II) have been computed using recent, consistent f-valucs (sec §6.1). 

b Indicatcs which lines of sight are used to compute the gas- phase Local Bubble (LBg) and gas-phase local disk (LDg) D/H values described in 
the text. INT indicates the intermediate regime, log N(H I) = 19.2-20.7. 

c D(Fe) and D(Si) refer to velocity Component C (Dupin & Cry 1998). N(H I) and N(D I) mostly in Component C. 

d -32 km s~ x component only. 

e D(Fe) and D(Si) refer to Component 1 (Howk et al. 1999), which is also called Component A (Shull & York 1977). 

f Abundances and depletions refer only to the H I component (Allen ct al. 1992). 

SN(H I) includes 2N(H 2 ) and (D/H) gas = [N(D I) + N(HD)]/[N(H I) + 2N(H 2 )]. 

h The Local Interstellar Cloud component only. For the BC component (D/H) gas = 51",! 1 ppm. 

References. — (1) Dring ct al. 1997. (2) Linsky ct al. 1995. (3) Wood ct al. 1996. (4) Wood ct al. 2000. (5) Piskunov ct al. 1997. (6) Wood ct al. 

2002. (7) Hcbrard et al. 2002. (8) Kruk ct al. 2002. (9) Lcmoine et al. 2002. (10) Vennes et al. 2000. (11) Lchncr ct al. 2002. (12) Oliveira et al. 

2003. (13) York & Rogerson 1976. (14) Sonncborn ct al. 2002. (15) Allen et al. 1992. (16) Cry ct al. 1985. (17) Friedman ct al. 2002. (18) Ferlet 



ct al. 1980. (19) York 1983. (20) Sonneborn et al. 2000. (21) Jenkins et al. 1999. (22) Laurent et al. 1979. (23) Wood et al. 2004. (24) Hoopes et 
al. 2003. (25) Willigcr ct al. 2005. (26) Fitzpatrick & Spitzer 1994. (27) Rcdficld & Linsky 2002. (28) Redfield & Linsky 2004a. (29) Dupin & Gry 
1998. (30) Lchncr ct al. 2003. (31) Bluhm et al. 1999. (32) van Steenberg & Shull 1988. (33) Mauche, Raymond, & Cordova 1988. (34) Spitzer 
& Fitzpatrick 1993. (35) Howk, Savage, & Fabian 1999. (36) Morton 1978. (37) Jenkins, Savage, & Spitzer 1986. (38) Prochaska et al. 2005. (39) 
Hebrard & Moos 2003. (40) Howk ct al. 2000. (41) Hebrard ct al. 2005. (42) Hebrard et al. 1999. (43) Friedman et al. 2006. (44) Olivcira ct al. 
2006. 



Table 3. COMPILATION OF D/H RATIOS AND METAL DEPLETIONS 



Target 


1 


b 


d 


logAf(HI) a (D/H) gas a 


(D/H) gas _ LB a 


D(Fc) a 


D(Si) a 




(dag) 


(deg) 


(PC) 


(ppm) 


(ppm) 







Sirius h 


227 


-09 


2.64 + 0.02 


17.60±g;i| 


16.2±?;* 




_1 11+0.12 
1 - il -0.14 


-0.63 + 0.16 


36 Oph 


358 


+07 


5.99 ± 0.04 


17.850 + 0.075 


15 ± 2.5 




-0.65 + 0.27 




e Eri 


227 


-48 


3.22 ± 0.01 


17.875 ± 0.035 


14.3 + 1.6 




-1.07 + 0.10 


-0.63 + 0.14 


31 Com 


115 


+89 


94 + 8 


17.884 + 0.03 


20.2 ± 1.9 




-0.78 + 0.10 




HZ 43 


54 


+84 


68 ± 13 


17.93 + 0.03 


16.6 ± 1.4 




-1.21 + 0.04 


-0.75 + 0.07 


e Ind 


336 


-48 


3.63 ± 0.01 


18.00 + 0.05 


16 + 2 








Procyon 


214 


+13 


3.50 ± 0.01 


18.06 + 0.05 


16 + 2 




-1.24 + 0.05 




P Cas 


118 


-03 


16.7 + 0.1 


18.132 + 0.025 


16.9 ± 1.6 




-1.22 ± 0.10 




HR 1099 


185 


-41 


29.0 ± 0.7 


18.131 ± 0.020 


14.6 ± 1.0 








G191-B2B 


156 


+07 


69 ± 15 


18.18 ± 0.09 


16.6 ± 4.1 




—0.58 ± 0.09 


^ 1 _i_0.17 
- L14 -0.16 


P CMa c 


226 


-14 


153 ± 15 


18.201°;^ 


> 16 




—0 fil + u - /u 


_i 1 1 ■>(! 

+0.04i°;20 


a Gem 


191 


+23 


37 ± 1 


18.201 ± 0.037 


13.8 ±2.1 








Capella 


163 


+05 


12.9 + 0.1 


18.239 + 0.035 


15.9 ± 1.5 




-1.20 + 0.04 


-0.75 + 0.06 


p bem 


i no 




in Q/i _l_ n nn 
lU.o4 ± U.Uy 


lo.zul it U.Uo / 


1/1 Q J- O O 

14. o ± Z.2 




i on J_ n nt; 
— i.zy ± U.Uo 




a in 




1 1 


iy. i ± u.o 


lo.oZ / it U.Uoo 


1 q q _l_ o n 
lo.o It Z.U 




i 1 q _|_ n on 
— ±. lo ± U.zU 




/3 Cet 


111 


— ol 


9Q A 4- fl 7 


1 Q QC. _|_ n nc; 
lO.OU HI U.UO 








tu.ou in Kj.^ti 


A And 


110 


— 15 


25.8 ± 0.5 


18.45 ± 0.075 


17 + 5 








Feige 24 


166 


—50 


74 ± 20 


18 47 + 03 


13 + 5 






—0 44 + 08 


WD 0621-376 


245 


-21 


78 + 23 


18.70 + 0.15 


14.1 + 6.0 




-1.21 + 0.20 




WD 2211-495 


346 


-53 


53 + 16 


18.76 + 0.15 


15.1 + 6.5 




-0.96 + 0.16 


-0.27 + 0.16 


WD 1634-573 


330 


-07 


37 + 3 


18.85 + 0.06 


15.8 + 2.5 




-1.22 + 0.09 


-0 60+ 14 


a Vir 


316 


+51 


80 + 6 


19.0 + 0.1 


15.8+^ 




-1 29+ - 10 

1 - zb '-0.12 


-0 S5+ 010 


GD 246 


87 


-45 


79 ± 24 


19.11 + 0.025 


15.1 ± 1.9 




-1.26 + 0.06 


-0 5o+°' 06 


A Sco d 


352 


-02 


216+42 


19.23 + 0.03 


7 6+1-8 

'■"-1.4 




-1.64 + 0.05 


-1.32 + 0.05 


P Cen 


312 


+01 


161 ± 15 


19.63 + 0.10 


11.7 + 7.5 


a - d -25.9 


-1.16 + 0.11 


-0.55 + 0.10 


7 2 Vel 


263 


-08 


258 ± 35 


19.710 + 0.026 


21.8 ± 2.1 


24.7 + 3.2 


-1.21 + 0.11 


-0.19 + 0.05 


a Cru 


300 


+00 


98 + 6 


19.85t° ;? 


12.6+j;? 


11 7+ 4 ' 7 
xl - ' -3.4 


-1 1 n +014 




BD+28°4211 


82 


-19 


104 ± 18 


19.846 + 0.018 


13.9 + 1.0 


13.4 ± 1.3 


-1.20 + 0.10 




Lan 23 


108 


-01 


122 ± 37 




21 9+ 4 ' 1 


no 

ZO.O_ 23 4 


i ti+0.07 
1 - ol -0.2S 




H Col c 


237 


-27 


4001^°° 


19.86 + 0.015 


6 9+ 6 ' 9 
°- a -1.8 




-1.18 + 0.02 


-0.27 + 0.02 


CPup 


256 


-05 


429 ± 94 


19.963 ± 0.026 


14.0 + 2.3 


13.7±l"i 


-l-2»_ .07 


-0 40+ 07 


TD1 32709 


233 


+28 


520 ± 90 


20.03 + 0.10 


18.6 ± 5.3 


19 1+ 6 ' 7 


-1.53 + 0.14 




WD 1034+001 


248 


+48 


155tf 3 


20.07 + 0.07 


21.4 + 5.3 


22.3 + 6.3 


-1.42 + 0.12 





Table 3 — Continued 



Target 


1 

(deg) 


b 

(deg) 


d 

(PC) 


logAf(H I) a 


(D/H) gas a 
(ppm) 


(D/H) gas _ LB a 
(ppm) 


D(Fe) a 


D(Si) a 


BD+39°3226 


65 


+29 


290+ 7u ° 


20 


.08 ± 0.09 


11.7 + 3.1 


11.2 + 3.4 


-1 


.38 ± 0.11 


-0.79 ± 0.21 


Feige 110 


74 


-59 


179±* B 


20 


i 4 +0.065 
■^-0.10 


21 4+ 5 - 7 


99 i+6.9 








7 Cas 


124 


-02 


188 ± 20 


20 


1fi +0.08 
1D -0.10 


9 8+ 2 - 7 


9-^2.5 


-1 


r-o+0.14 
oo -0.18 




i Ori 


210 


-20 


407 ± 127 


20 


i c+0.06 
10 -0.07 


14.1 ± 2.8 


13.9 + 3.1 


-1 


■^ 1 -0.16 


-0.51 ± 0.07 


<5 Ori 


204 


-18 


281 ± 65 


20 


19 ± 0.03 


7 4+1-2 


D '°-1.0 


-1 


.56 ± 0.04 




Car f 


290 


-05 


135 ± 9 


20 


.28 ± 0.10 


5 n+ 2 - 9 




-1 


.54 + 0.11 


_, 41 +0.26 

i -^ i -0.12 


e Ori 


205 


-17 


412 ± 154 


20 


4E+O.O8 




5.8±i; 8 5 


-1 


70 +0.14 
' u -0.12 


-0 94+°' 12 
u - a4 -0.14 


PG 0038+199 g 


120 


-43 


„ q7 +l64 
za '-104 


20 


.48 ± 0.04 


19.1 ± 2.6 


19.2 + 2.8 


-1 


ci +0.05 
oi -0.04 




LSE 44 


313 


+13 


554 ± 66 


20 


cn+0.10 
oz -0.18 




99 7+12.8 
""■'-6.6 






JL 9 


323 


-27 


590 ± 160 


20 


.78 ± 0.05 


10.0 + 1.9 


9.8 ± 2.0 


-1 


.54 + 0.10 




HD 195965 


86 


+05 


794 ± 200 


20 


.95 ± 0.025 


8.5 + 1.6 


8.4 + 1.6 


-1 


.59 ± 0.03 




LSS 1274 


277 


-05 


580 ± 100 


20 


.98 ± 0.04 


7.6 + 1.9 


7.5 ± 1.9 


-1 


.62 ± 0.08 




HD 191877 


62 


-06 


2200 ± 550 


21 


.05 ± 0.05 


7 S+ 2 - 4 
'• 8 -1.5 


7.6±ll 


-1 


.55 ± 0.05 




HD 90087 s 


285 


-02 


2740 ± 800 


21 


.22 ± 0.05 


8.7 + 1.7 


8.6 ± 1.7 


-1 


,45 ± 0.05 





a Quoted uncertainties assumed to be lcr (see text). N(H I) is the most recent published value. N(D I) is not listed when not 
in original paper. D/H, D(Fe), and D(Si) arc computed using the tablulated value of N(H I). 

b Indicates which lines of sight arc used to compute the gas-phase Local Bubble (LBg) and gas-phase local disk (LDg) 
D/H values described in the text. INT indicates the intermediate regime, log N(H I) = 19.2-20.7. 

c D(Fe) and D(Si) refer to velocity Component C (Dupin & Gry 1998). N(H I) and N(D I) mostly in Component C. 

d — 32 km s^ 1 component only. 

e D(Fe) and D(Si) refer to Component 1 (Howk ct al. 1999), which is also called Component A (Shull & York 1977). 
Abundances and depletions refer only to the H I component (Allen ct al. 1992). 
S N(H I) includes 2N(H 2 ) and (D/H) gas = [N(D I) + N(HD)]/[N(H I) + 2N(H 2 )]. 
h The Local Interstellar Cloud component only. For the BC component (D/H) gas = ppm. 
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Table 4. ESTIMATES OF (D/H) gas FROM (D/0) gas x (O/H) 



Reference 


Oxygen 


AT 1 

IN umber 


Oxygen LOS 


( r\ /tt\ 
(0/H)gas 


(D/H) gas a 




Data oource 


OI LUb 


distance (pc) 


(ppm) 


(ppm) 


HebrarH Rr Moos (2003^1 












cite Meyer (2001) 


GHRS, STIS 


29 


150-5010 


343 ± 15 


13.2 ±0.8 


Andre et al. (2003) 


FUSE 


19 


840-5010 


408 ± 13 


15.7 ±0.9 


include Meyer et al. (1998) 


GHRS 


32 


150-5010 


377 ± 10 


14.5 ±0.7 


include 5 dense LOS 


GHRS 


37 


150-5010 


362 ± 11 


13.9 ±0.7 


Oliveira et al. (2005) 


FUSE 


11 


2.6-128 


345 ± 19 


13.2 ±0.9 


Frisch & Slavin (2003) b 


Models 2,8 




3 


380 ± 30 


14.6 ± 1.3 


Asplund et al. (2004) 


solar c 






458 ± 53 


17.6 ±2.1 


if 120 ppm in dust 








338 ± 53 


13.0 ±2.0 


Melendez (2004) 


solar c 






390 ± 63 


15.0 ±2.4 


Wood et al. (2004) 






3.2-94 




15.6 ±0.4 



a Assuming (D/0) gas = (3.84 ±0.16) x 10~ 2 (la) (Hebrard & Moos 2003). 

b Frisch & Slavin (2003) do not assign an error to (0/H) gas . We arbitrarily assign an error of 
8% based on the discussion in their paper. 

c Column 5 lists (O/H) , the total abundance ratio in the solar photosphere. 



Table 5. CORRELATION PARAMETERS 



Correlation 3 


log AT(H I) 


N 


a 


±1<T 


b 


±1(7 


Spearman 


Significance 


(D/H) gas vs. D(Fc) 


all 


38 


30.00 


1.60 


13.04 


1.15 


0.483 


0.0021 (2.9cr) 


(D/H) gas vs. D(Fe) 


HST,FUSE,IMAPS 


29 


29.14 


1.73 


12.24 


1.30 


0.398 


0.0324 (2.1a) 


(D/H) gas vs. D(Fe) 


> 19.0 


24 


35.04 


3.54 


16.68 


2.37 


0.534 


0.0072 (2.6cr) 


(D/H) gas vs. D(Si) 


all 


20 


20.09 


0.81 


9.72 


1.05 


0.325 


0.157 (1.4cr) 


(D/H) gas vs. D(Si) 


> 19.0 


11 


19.27 


0.86 


10.79 


0.984 


0.664 


0.026 (2.1cr) 


(D/Fe) gas vs logiV(H I) 


all 


38 


-0.59 


0.12 


0.081 


0.006 


0.516 


0.0009 (3.1a) 


(D/Fe) gas vs log N(R I) 


> 19.2 


22 


1.31 


0.47 


-0.011 


0.023 


0.0936 


0.679 (0.4cr) 


(D/H) gas _ LB vs. logT 01 


all 


16 


-10.66 


4.45 


10.54 


2.22 


0.564 


0.023 (2.2cr) 



Correlation equations are of the form y = a + bx, where a is the intercept and b is the slope. 
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Table 6. MOLECULAR HYDROGEN TEMPERATURE MEASUREMENTS 



Target 


logiV(H I) a 


logAf(H 2 ) a 


log/(H 2 ) a 


(D/H) gas -LB a 
(ppm) 


T i a 
(K) 


Refs. 


7 2 Vcl d 


19.710 ±0.026 


14.23 


-5.18 


24.7 ±3.2 


c-t-o + llSO 
ddo -292 


1,2 


BD +28°4211 


19.846 ±0.018 


K 1O+0.20 
1J ' lo -0 .08 


_4 4-3+0.20 


13.4 ± 1.3 




11 


Lan 23 


19-89±g:gS 


15.11 


-4.48 


no 


286±Sg° 


4 


H Col d 


19.86 ±0.015 


15.51 


-4.09 




98111 


2,5 


C Pup d 


19.963 ± 0.026 


14.45 


-5.21 


13.7±|-i 


382tf 5 ° 


1,2 


TD1 32709 c 


20.03 ±0.10 


14.481°;!? 


r ncr+0.15 
°- zo -0.14 


i9-i±l:I 


292 ± 83 


14 


WD 1034+001 c 


20.07 ±0.07 


15.721°;" 


-4 05+ 14 


22.3 ±6.3 


341 ± 75 


14 


BD+39°3226 


20.08 ± 0.09 


IO. Od_Q q 7 


-4.13 ±0.11 


11.2 ±3.4 


104 ± 27 


14 


i Ori d 


20 16+ 06 


14.69 


-5.17 


13.6+J; 


481+gl 


2,6 


8 Ori 


20.19 ±0.03 


14.74 


-5.16 


6-5±i;l 


598 -243 


7,13 


6 Car 


20.28 ±0.10 


15.02 


-4.96 


4.0+1° 


96 ±9 


10 


e Ori 


20.451°;?* 


16.28 


-3.89 


5 8+ 1 - 8 


89l|° 


6,13 


PG 0038+199 b 


20.48 ± 0.04 


19.33 ±0.02 


-0.85 ±0.05 


19.2 ± 2.8 


139 ±9 


8 


JL 9 b 


20.81 ±0.05 


19.25 ±0.015 


-1.26 ±0.06 


9.3 ±2.0 


89 ± 3 


9 


LSS 1274 b 


20.99 ±0.04 


19.10 ±0.02 


-1.59 ±0.05 


7.3 ± 1.9 


64 ±2 


9 


HD 90087 b 


21.22 ±0.05 


19.92 ±0.02 


-1.00 ±0.05 


8.6 ±1.7 


79 ± 2 


12 



a Quoted uncertainties are la. /(H 2 ) = 2Af(H 2 )/[2A r (H 2 ) + N(U I)]. 
T i = 171//n[9Af(J = 0)/AT(J = 1)] K. 

b 7V(H I) includes 2W(H 2 ) and (D/H) gas _ LB = [JV(D I) + JV(HD)]/[JV(H I) ± 2JV(H 2 )]. 

c Tq2 is given instead of Toi which is indeterminate. 

d Errors in N(J = 0) and N(J = 1) arc assumed to be 0.08 (la). 

References. — (1) Sonneborn ct al. 2000. (2) Savage et al. 1977. (3) Bluhm ct al. 1999. (4) Oliveira 
ct al. 2003. (5) York & Rogerson 1976. (6) Laurent ct al. 1979. (7) Jenkins et al. 1999. (8) Williger et al. 
2005. (9) Wood ct al. 2004. (10) Allen et al. 1992. (11) Sonneborn ct al. 2002. (12) Hcbrard ct al. 2005. 
(13) Jenkins et al. 2000b. (14) Oliveira ct al. 2006. 
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Fig. 1. — Plot of the gas-phase D/H values (not corrected for the Local Bubble foreground) 
vs. the hydrogen column density for the 47 lines of sight listed in Table 3. The symbols 
for each data point indicate the spacecraft that observed the line of sight. Error bars are 
±1 a. Vertical dashed lines indicate lines of sight inside the Local Bubble [logiV(H I) < 19.2 
cm -2 ] and the intermediate region [logiV(H I) between 19.2 and 20.7]. The solid horizontal 
lines indicate the mean values of (D/H) gas for the low and high iV(H I) regions, and the 
dotted horizontal lines indicate the ±1 a about the mean value. The dash-dot horizontal 
line is the mean value of (D/H) gas _LB for the highest five points in the intermediate region 
(7 2 Vel, Lan 23, WD 1034+001, Feige 110, and LSE 44) after subtracting the foreground 
Local Bubble contributions to the hydrogen and deuterium column densities along the lines 
of sight. See discussion in §7. 
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Fig. 2. — (D/H) gas vs. (0/H) gas with la errors (see text). The data do not show anticorrelation of (D/H) 
with (O/H), as would result from variable astration or infall, such as model C-I from Chiappini et al. (2002). 
The full range of D/H is seen for sightlines with (0/H) gas between 300 and 390 ppm. Depletion of D into 
interstellar dust appears to be the most viable explanation for the observed variations in (D/H) gas . 
Sightlines: WD 1634-573 (Wood et al. 2002); CD 246 (Olivcira et al. 2003); BD+284211 (Sonncborn et al. 
2002; Spitzer et al. 1974); Fcigc 110 (Fricdmann ct al. 2002; Hebrard et al. 2005); 7 Cas (Ferlet et al. 1980; 
Meyer et al. 1988; Meyer 2001); S Ori (Jenkins et al. 1999; Meyer et al. 1998; Meyer 2001); 1 Ori (Laurent et 
al. 1979; Meyer et al. 1998; Meyer 2001); HD 195965 and HD 191877 (Hoopes et al. 2003); HZ 43A: (Kruk 
et al. 2002); G 191-B2B (Lemoine et al. 2002); LSS 1274 (Wood et al. 2004); Hebrard & Moos (2003). 
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Fig. 3. — Plot of gas-phase D/H values toward stars vs. the depletion of iron for the 38 lines 
of sight listed in Table 3 (ignoring upper limits). The symbols are the same as in Figure 1, 
and the error bars are ±1 a. The solid line is the least-squares linear fit to the data weighted 
by the inverse errors. See §6 for a discussion of the error analysis technique. Five lines of 
sight with (D/H) gas well above the linear fit are identified. The dashed line has the same 
slope as the solid line but is displaced upward by 8.5 ppm. See discussion in §6.3.4. 
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Fig. 4. — Plot of gas-phase D/H values toward stars vs. the depletion of silicon for the 
20 lines of sight (ignoring upper limits) listed in Table 3. The symbols are the same as in 
Figure 1, and the error bars are ±1 a. The solid line is the least-squares linear fit to the 
data (except for upper limits) weighted by the inverse errors. Lines of sight to three stars 
are identified. 
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Fig. 5. — Plot of gas-phase D/H values toward stars (after correction for the Local Bubble 
foreground) vs. T 01 , the excitation temperature of molecular hydrogen between the J=0 and 
J=l rotational levels. The solid line is the least-squares linear fit to the 16 data points (see 
Table 6) weighted by the inverse errors. Two outliers are identified. 
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Fig. 6. — Plot of gas-phase D/H values toward stars vs. the depletion of iron for the 24 lines 
of sight listed in Table 3 with logiV(H I) > 19.0. The symbols are the same as in Figure 1, 
and the error bars are ±1 a. The solid line is the least-squares linear fit to the data weighted 
by the inverse errors. The dashed line is the fit to the 38 lines of sight covering the full range 
of log iV(H I) as shown in Figure 3. 
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Fig. 7. — Plot of gas-phase D/H values toward stars vs. the depletion of silicon for the 
11 lines of sight listed in Table 3 with logiV(H I) > 19.0. The symbols are the same as in 
Figure 1, and the error bars are ±1 a. The solid line is the least-squares linear fit to the data 
weighted by the inverse errors. The dashed line is the fit to the 20 lines of sight covering the 
full range of logiV(H I) as shown in Figure 4. 
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Fig. 9. — Column density ratios log[iV(D I)/iV(Fe II)] vs. logiV(H I) for all lines of sight 
with symbol coding the same as in Figure 1. The solid line is a weighted least-squares fit to 
all of the data points, and the dashed line is a similar fit to only the data beyond the Local 
Bubble (logiV(H I) > 19.2). 



